NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


PASSIVE  DETECTION  OF  GASES  IN  THE  ATMOSPHERE. 
CASE  STUDY:  REMOTE  SENSING  OF  SO,  IN  THE  UV 
USING  LINUS 

by 

Anastasios  G.  Halvatzis 

Thesis  Advisor: 

Second  Reader: 

December  2002 

Richard  C.  Olsen 

Richard  M.  Harkins 

Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


_ REPORT  DOCUMENTATION  PAGE _ |  Form  Approved  OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including 
the  time  for  reviewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and 
completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any 
other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington 
headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project 
(0704-0188)  Washington  DC  20503. _ 


1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

December  2002 

3.  REPORT  TYPE  AND  DATES  COVERED 

Master’s  Thesis 

4.  TITLE  AND  SUBTITLE:  Passive  Detection  of  Gases  in  the  Atmosphere.  Case 
Study:  Remote  Sensing  of  SO?  in  the  UV  Using  LINUS 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S):  Anastasios  G.  Halvatzis 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

N/A 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES:  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
jrolic^o^jrositioir^ofthe^ejrartmeinofDefense^oMlm^S^^Goveniment^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited 


13.  ABSTRACT  (maximum  200  words) 


An  imaging  UV  spectrometer  was  used  to  study  sulfuric  plumes  at  Lassen  Volcanic  National  Park,  in  an 
effort  to  identify  and  quantify  S02  presence.  The  NPS  instrument,  LINUS  (Lineate  Imaging  Near -Ultraviolet 
Spectrometer)  was  taken  to  Lassen  on  September  13  and  14,  2002.  Data  taken  there  are  compared  to  laboratory 
measurements  of  S02  subsequently  run  at  NPS.  These  data,  along  with  additional  measurements  of  platinum 
discharges  for  wavelength  calibration,  allow  for  comparison  with  atmospheric  modeling  calculations. 
Observations  were  modeled  with  the  standard  MODTRAN  code.  Comparisons  between  simulated  and  measured 
data  showed  minor  indications  of  S02  in  the  Lassen  data.  The  gas  concentration  was  estimated  to  be  less  than 
lOppmv. 


14.  SUBJECT  TERMS :  Gas  Detection,  Sulfur  Dioxide,  Remote  Sensing,  Ultraviolet  (UV),  LINUS, 
Spectral  Imager,  MODTRAN 

15.  NUMBER  OF 
PAGES 

135 

16.  PRICE  CODE 

17.  SECURITY 

18.  SECURITY 

19.  SECURITY 

20.  LIMITATION 

CLASSIFICATION  OF 

CLASSIFICATION  OF  THIS 

CLASSIFICATION  OF 

OF  ABSTRACT 

REPORT 

PAGE 

ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

UL 

NSN  75404)1-280-5500 


Standard  Form  298  (Rev.  2  -89) 
Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


PASSIVE  DETECTION  OF  GASES  IN  THE  ATMOSPHERE. 
CASE  STUDY:  REMOTE  SENSING  OF  S02  IN  THE  UV  USING  LINUS 

Anastasios  G.  Halvatzis 
Major,  Hellenic  Air  Force 
B.S.,  Hellenic  Air  Force  Academy,  1988 
B.S.,  University  of  Macedonia  in  Thessaloniki,  1995 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  PHYSICS 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  2002 


Author:  Anastasios  G.  Halvatzis 


Approved  by:  Richard  C.  Olsen 

Thesis  Advisor 


Richard  M.  Harkins 
Second  Reader 


William  B.  Maier  II 
Chairman,  Department  of  Physics 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


ABSTRACT 


An  imaging  UV  spectrometer  was  used  to  study  sulfuric  plumes  at  Lassen 
Volcanic  National  Park,  in  an  effort  to  identify  and  quantify  S02  presence.  The  NPS 
instrument,  LINUS  (Lineate  Imaging  Near -Ultraviolet  Spectrometer),  was  taken  to 
Lassen  on  September  13  and  14,  2002.  Data  taken  there  are  compared  to  laboratory 
measurements  of  S02  subsequently  run  at  NPS.  These  data,  along  with  additional 
measurements  of  platinum  discharges  for  wavelength  calibration,  allow  for  comparison 
with  atmospheric  modeling  calculations.  Observations  were  modeled  with  the  standard 
MODTRAN  code.  Comparisons  between  simulated  and  measured  data  showed  minor 
indications  of  S02  in  the  Lassen  data.  The  gas  concentration  was  estimated  to  be  less  than 
lOppmv. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

H.  DEFINITIONS  AND  BASIC  PHYSICS . 3 

A.  REMOTE  SENSING . 3 

B.  ELECTROMAGNETIC  RADIATION . 4 

1.  Electromagnetic  Spectrum . 4 

2.  Blackbody  Radiation . 5 

3.  Interactions  with  the  Atmosphere . 7 

a.  Scattering. . 9 

b.  Transmission . 9 

c.  Absorption . 9 

d.  Reflection . 10 

III.  SULFUR  DIOXIDE . 11 

A.  GENERAL . 11 

B.  ABSORPTION  PROCESSES  IN  THE  VISIBLE  AND  UV . 11 

C.  BEER’S  LAW  -  ABSORPTION  CROSS  SECTION . 12 

D.  S02  ABSORPTION  CROSS  SECTIONS  IN  THE  NEAR  UV . 13 

IV.  UV  REMOTE  SENSING  OF  VOLCANIC  SO  2 . 19 

A.  GENERAL . 19 

B.  TOMS . 19 

C.  BREWER  SPECTROMETER . 22 

D.  COSPEC . 24 

1.  Brief  Description  of  Instrument  and  Its  Operation . 24 

2.  Calculations  and  Limitations  of  the  SO  2  Flux  Determination . 26 

3.  S02  Flux  Measurement  Techniques  Using  COSPEC . 29 

E  FLYSPEC . 31 

F.  SUMMARY . 33 

V.  LINEATE  IMAGING  NEAR  ULTRA  VIOLET  SPECTROMETER 

(LINUS) . 37 

A.  GENERAL . 37 

B.  DEVICE  DESCRIPTION  AND  OPERATION . 37 

1.  Hardware . 37 

a.  Host  Computer. . 38 

b.  Scanning  system . 38 

c.  Camera  system . 38 

d.  Optics . 39 

2.  Software . 40 

3.  LINUS  Operation . 40 

VI.  EXPERIMENTAL  DATA  COLLECTION . 43 

A.  GENERAL . 43 

vii 


B.  BRIEF  VOLCANIC  HISTO  RY  OF  LASSEN  PARK . 44 

C.  DATA  ACQUISITION . 45 

1.  First  Trip . 45 

2.  Second  Trip . 45 

D.  RAW  DATA . 49 

VII.  MODELING  WITH  MODTRAN . 53 

A.  GENERAL . 53 

B.  MODTRAN4 . 53 

C.  SELECTION  OF  MODTRAN4  INPUT  CARDS  AND  VARIABLE 

VALUES . 55 

1 .  Selection  of  Input  Cards . 56 

2.  Selection  of  Variable  Values . 56 

a.  Card  1. . 57 

b.  CardlA . 59 

c.  Card  2 . 61 

d.  Card  2C . 61 

e.  Cards  2C1  and  2C2 . 61 

f.  Card  2 . 63 

g.  Card  3A1 . 63 

h.  Card  3A2. . 64 

i.  Card  4. . 64 

D.  SIMULATION  NOTES . 64 

1.  Variable  Values  Verification . 64 

2.  Computational  Time . 65 

3.  Spectral  Resolution . 65 

VIII.  DATA  ANALYSIS . 67 

A.  SIMULATED  DATA . 67 

1.  LINUS  UV  Filter . 67 

2.  S02  Presence  Effects  and  Spectral  Signature . 68 

a.  Figures  8.2  through  8.4. . 68 

b.  Figures  8.5  and  8.6 . 71 

c.  Figures  8.7  and  8.8 . 74 

3.  Spectral  Profile  Structure  Investigation . 77 

B.  EXPERIMENTAL  DATA . 79 

1.  Wavelength  and  S02  Calibration  of  LINUS . 79 

2.  Data  Analysis . 82 

a.  Individual  Image  Analysis . 82 

b.  Ensemble  Analysis . 87 

C.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  SIMULATED 

DATA . 88 

D.  SUMMARY . 93 

IX.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS  . 97 

APPENDIX  A.  MODTRAN  INPUT  CARDS  AND  VARIABLES . 99 

viii 


APPENDIX  B.  EXAMPLES  OF  MODTRAN  INPUT  CARDS . 101 

APPENDIX  C.  SELECTED  SPECTRAL  PROFILES  FROM 

EXPERIMENTAL  DATA  . 107 

LIST  OF  REFERENCES . Ill 

INITIAL  DISTRIBUTION  LIST . 115 


IX 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


x 


LIST  OF  FIGURES 


Figure  2.1.  Spectral  imager  classifications.  (From:  Multispectral  Users  Guide,  1995) . 4 

Figure  2.2.  Electromagnetic  spectrum.  (From:  Barrett  &  Curtis,  1982) . 5 

Figure  2.3.  Blackbody  radiation  curves  for  various  temperatures.  As  shown  by  (a)  and 

(b)  as  the  temperature  increases  there  is  a  significant  change  of  scale  and 
the  maximum  irradiance  occurs  at  shorter  wavelengths.  (From:  Elachi, 

1987) . 6 

Figure  2.4.  Atmospheric  layers.  (From:  Lenoble,  1993) . 8 

Figure  2.5.  Interaction  of  EM  radiation  with  the  atmosphere.  (From:  Gibson,  2000) . 8 

Figure  2.6.  Solar  irradiance  approximated  by  a  blackbody,  at  the  top  of  the 

atmosphere  (TOA)  and  at  sea  level.  (From:  Elachi,  1987) . 10 

Figure  3.1.  Absorption  spectrum  of  SO2.  (From:  Atkins,  1982) . 12 

Figure  3.2.  SO2  absorption  cross  section  from  290  to  317nm.  (From:  Brassington, 

1981) . 15 

Figure  3.3.  Comparison  of  measurements  between  Brassington  (solid  line)  and 

Thomson  et  al.  (dotted  line).  (From:  Brassington,  1981) . 16 

Figure  3.4.  Comparison  of  measurements  between  Brassington,  Thomson  et  al.  and 

Woods  et  al.  for  the  300.03nm  absorption  peak.  (From:  Brassington,  1981)  .  16 
Figure  3.5.  High  precision  SO2  absorption  cross-section  peaks  at  298  and  300nm  at 


20"C.  The  dotted  line  is  with  l.OlxlO5  Pa  of  nitrogen  buffer  gas.  The  solid 
lines  are  is  for  SO2  alone.  The  original  article  refers  to  the  solid  line  as  the 
one  with  the  buffer  gas  and  the  dotted  as  SO2  alone,  but  this  is  obviously 


opposite  since  the  buffer  gas  smoothes  out  the  spectral  details.  (From: 

Woods  et  al.,  1980) . 17 

Figure  3.6.  High  resolution  absorption  cross  section  of  SO 2  at  20°C  with  nitrogen 
buffer  gas  over  the  spectral  region  297.8-300.3nm.  (From:  Woods  et  al., 

1980) . 17 

Figure  3.7.  SO2  absorption  cross  section  for  (a)  295  K  and  (b)  210  K.  (From:  Thomas 

etal.,  1987) . 18 

Figure  4.1.  Absorption  cross  sections  of  atmospheric  UV  absorbers.  (From: 

Strickland,  1999) . 20 

Figure  4.2.  False-color  image  of  S02  from  the  June  16,  1991  eruption  of  Mt. 

Pinatubo,  Philippines,  produced  from  TOMS  data.  The  triangular  symbol 
marks  the  location  of  the  source  volcano.  (From: 

<http://skye.gsfc.nasa.gov>,  2002) . 21 

Figure  4.3.  Brewer  Station  on  the  roof  of  the  Environment  Canada  Building,  Toronto. 

(From:  <http://www.ns.ec.gc.ca>,  2002) . 23 

Figure  4.4.  Example  of  Brewer  station  output.  (From:  WOUDC,  2002) . 23 

Figure  4.5.  Close-up  view  of  COSPEC.  (From:  USGS  website,  2002) . 24 

Figure  4.6.  Internal  design  of  COSPEC.  (From:  COSPEC  IV  Manual,  1976) . 25 

Figure  4.7  Example  of  COSPEC  output  to  a  paper  chart  recorder.  (From: 

<http://www.gsj.jp>,  2002) . 26 


Figure  4.8.  Example  of  long-term  COSPEC  calculation  results  in  metric  tons/day 
(mt/d)  of  SO?  Kilauea  Summit,  Hawaii  SO?  emission  rates.  (From:  Elias  et 

al.,  1998).....". . 27 

Figure  4.9.  Example  of  estimated  uncertainties  of  SO2  flux  calculations  based  on 
vehicle -mounted  COSPEC.  The  black  vertical  bars  represent  the  standard 
deviation  of  all  traverses  on  a  single  day.  Kilauea’ s  East  Rift  Zone, 

Hawaii.  (From:  Elias  et  al.,  1998) . 28 

Figure  4.10.  COSPEC  mounted  on  vehicle,  (From:  <http://www.volcanoes.ca>,  2002)  ...30 

Figure  4. 1 1.  COSPEC  on  airplane  and  on  ground.  (From:  USGS  website,  2002)  . 30 

Figure  4.12.  Concurrent  data  collection  with  COSPEC  and  FLYSPEC.  (From:  Mares, 
2002)  32 

Figure  4.13.  FLYSPEC  sample  data  frame  collected  at  Kilauea  Volcano,  Hawaii  in 

March  2002.  (From:  Mares,  2002) . 32 

Figure  4.14.  Comparison  between  COSPEC  and  FLYSPEC  data  taken  concurrently  in 

March  2002  at  Kilauea  Volcano,  Hawaii.  (From:  Mares,  2002) . 33 

Figure  4.15.  UV  sensors  used  in  the  detection  of  volcanic  S02-  (From:  Krueger  at  al., 
2000)  34 

Figure  4.16.  Absorption  spectra  of  SO2  and  O3  and  wavelength  bands  used  by  the  UV 

detectors.  (From:  Krueger  at  al.,  2000) . 35 

Figure  5.1.  LINUS  architecture.  (From:  Kompatzki,  2000) . 37 

Figure  5.2.  LINUS  optical  layout.  (From:  Kompatzki,  2000) . 39 

Figure  5.3.  Spectral  response  for  the  300nm  UV  filter.  (From:  Filter  Specs) . 39 

Figure  5.4.  Pictures  of  LINUS.  On  the  left  side  is  the  electronics  section  with  the 

computer  the  amplifier  and  the  camera  pulser  and  cooler.  On  the  right  side 
is  the  optics  section,  which  also  includes  the  camera  body  and  the  scanner 

miiTor  assembly . 40 

Figure  5.5.  LINUS  spectral  imaging  schematic  diagram . 42 

Figure  5.6.  Data  set  produced  by  LINUS  output . 42 

Figure  6.1.  Map  of  Lassen  Volcanic  National  Park.  The  square  on  the  right  is 

magnification  of  the  area  included  in  the  left  box.  Locations  1,  2  and  3  are 
the  spots  where  the  measurements  took  place.  (After:  USGS  website, 

2002) . 43 

Figure  6.2.  Some  of  the  Sulfur  Works  vents . 44 

Figure  6.3.  The  May  22,  1915  explosive  eruption  of  Lassen  Peak.  (From:  USGS 

website,  2002) . 45 

Figure  6.4.  Sulfur  Works  as  seen  from  location  1 . 47 

Figure  6.5.  LINUS  setup  at  location  2 . 48 

Figure  6.6.  Screenshot  of  LINUS  monitor  at  the  time  of  data  file  acquisition.  This  is 

how  the  raw  data  looked  like  to  the  operators . 50 

Figure  6.7.  The  previously  presented  data  file  in  pseudocolor,  which  gave  to  Linus 
operators  a  first  feeling  about  the  intensity  variations  along  the  spectrum 

contained  in  the  image . 51 

Figure  6.8.  A  first  indication  of  the  spectral  profile  of  the  recorded  image  by 
averaging  the  middle  256  lines.  Vertical  axis  units  are  Digital  Numbers 
(DN) . 52 

xii 


Figure  7.1.  MODTRAN  results  for  single  and  multiple  scattering . 58 

Figure  7.2.  Ratio  of  single  divided  by  multiple  scattering  from  MODTRAN  results . 58 

Figure  7.3.  Effect  of  using  the  DISORT  algorithm  in  the  spectral  region  of  interest . 59 

Figure  7.4.  Solar  irradiance  comparison  for  5cm  1  (LSUN=F)  and  1cm  1  (LSUN=T) 

spectral  resolution.  Results  for  Lassen  Park  with  MODTRAN  standard 

atmosphere  2 . 60 

Figure  7.5.  SO2  concentration  as  a  function  of  altitude  for  the  MODTRAN  standard 

atmosphere  2 . 63 

Figure  8.1.  The  continuous  line  shows  the  Gaussian  fit  of  the  LINUS  filter  response. 

Diamonds  show  data  points  taken  from  the  real  filter  response  curve 

presented  in  figure  5.3 . 67 

Figure  8.2.  Simulated  sky  radiance  spectral  profiles  for  various  SO2  concentrations  at 

Sulfur  Works  without  LINUS  filter  effect.  Resolution  0.08nm . 69 

Ligure  8.3.  Simulated  sky  radiance  spectral  profiles  for  various  SO2  concentrations  at 
Sulfur  Works,  without  LINUS  filter  effect,  at  smaller  scale  than  figure  8.2. 

Resolution  0.08nm . 70 

Figure  8.4.  Simulated  sky  radiance  spectral  profiles  for  various  SO2  concentrations  at 

Sulfur  Works  without  LINUS  filter  effect.  Resolution  0.5nm . 70 

Figure  8.5.  Simulated  spectral  profiles  for  various  SO2  cone  entrations  at  Sulfur  Works 

with  LINUS  filter  effect.  Resolution  0.08nm . 72 

Figure  8.6.  Simulated  spectral  profiles  for  various  SO2  concentrations  at  Sulfur  Works 

with  LINUS  filter  effect.  Resolution  0.5nm . 73 

Ligure  8.7.  Ratios  of  SO2  concentrations  to  regular  atmosphere  at  0.08nm  resolution . 75 

Ligure  8.8.  Ratios  of  S02  concentrations  to  regular  atmosphere  at  0.5nm  resolution . 76 

Ligure  8.9.  Spectrum  of  regular  atmosphere  and  atmosphere  with  no  03 . 77 

Ligure  8.10.  Ratio  of  radiances  in  log-scale  of  regular  atmosphere  to  atmosphere  with 

no  O3 . 78 

Ligure  8.11.  Linear  ratio  of  radiances  of  regular  atmosphere  to  atmosphere  with  no  0>...  78 

Ligure  8.12.  Spectrums  of  solar  irradiance  out  of  the  atmosphere  and  regular 

atmosphere  radiance  at  2.13km  altitude . 79 

Ligure  8.13.  Results  of  wavelength  calibration  of  LINUS  using  a  platinum  lamp . 80 

Ligure  8. 14.  S02  calibration  results  of  LINUS . 81 

Ligure  8.15.  Plot  of  the  averages  from  figures  C.2  through  C.5 . 83 

Ligure  8.16.  Comparison  of  spectral  profiles  with  strong  structure  over  the  Sulfur 
Works  vents  (possible  SO2  presence)  with  Lassen  Peak  (clear 

atmosphere) . 84 

Ligure  8.17.  Ratios  of  strong  Sulfur  Works  spectra  (possible  SO  2  existence)  with  strong 

Lassen  Peak  spectra  (clear  atmosphere) . 85 

Ligure  8.18.  Ratios  of  the  Sulfur  Works  averages  (blue,  green)  to  the  Lassen  Peak 

average  (red)  shown  in  figure  8.15  above . 86 

Ligure  8. 19.  Ratios  of  Sulfur  Works  (possible  SO2  existence)  with  Lassen  Peak  spectra 
(clear  atmosphere).  Spectral  profiles  that  do  not  have  strong  structure  were 

used . 86 

Ligure  8.20.  Summary  of  the  experimental  measurements  with  a  sim  ulated  clear 

atmosphere  and  SO2  calibration  spectral  profiles . 87 


xiii 


Figure  8.21.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere.  The 

simulated  data  are  in  0.08nm  resolution . 89 

Figure  8.22.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere.  The 

simulated  data  are  in  0.5nm  resolution . 89 

Figure  8.23.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere  and 
simulated  with  lOppmv  SCF  atmosphere.  The  simulated  data  are  in  0.5nm 

resolution . 90 

Figure  8.24.  Ratios  of  data  taken  over  the  Sulfur  Works  vents  (SW_080  and  SW_206) 
and  clear  atmosphere  at  Lassen  Peak  (SW_503)  divided  by  simulated  clear 

atmosphere . 91 

Figure  8.25.  Inverted  ratios  of  simulated  clear  atmosphere  and  lOppmv  of  SCL- 
containing  atmosphere  divided  by  clear  atmosphere  at  Lassen  Peak 

(SW_503) . 92 

Figure  8.26.  Ratios  of  selected  experimental  to  simulated  data . 92 

Figure  8.27.  Comparison  between  MODTRAN  and  SCL  calibration  absorption 

characteristics . 93 

Figure  C.l.  Plot  of  the  spectral  profiles  of  all  selected  images . 107 

Figure  C.2.  Plot  of  spectral  profiles  of  one  group  of  selected  images  at  Sulfur  Works...  108 

Figure  C.3.  Plot  of  spectral  profiles  of  another  group  of  selected  images  at  Sulfur 

Works . 108 

Figure  C.4.  Plot  of  spectral  profiles  of  selected  images  at  Lassen  Peak . 109 

Figure  C.5.  Plot  of  spectral  profiles  of  measurements  taken  with  the  aperture  of 

LINUS  covered . 109 


xiv 


LIST  OF  TABLES 


Table  2.1.  Atmospheric  constituents.  (From:  Durkee,  1999) . 7 

Table  4.1.  Satellite  missions  with  TOMS  instruments.  (From:  Krueger  et  al.,  2000) . 22 

Table  4.2.  Calculation  of  S02  flux  from  COSPEC  output.  (From:  Stoiber  et  al.,  1983). .27 

Table  4.3.  Limitations  of  the  COSPEC  calculations.  (From:  Stoiber  et  al.,  1983) . 28 

Table  7.1.  Molecular  species  for  user  defined  atmospheric  profiles.  S02  is  number  9. 

(After  MODTRAN4  User’s  Manual,  1999) . 55 


xv 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xvi 


ACKNOWLEDGMENTS 


I  would  like  to  express  my  acknowledgments  to  the  people  that  directly  or 
indirectly  contributed  to  the  creation  of  this  thesis.  First  of  all,  many  thanks  to  Professor 
R.  C.  Olsen  for  his  help  and  guidance.  His  skills  and  experience  with  IDL  were 
invaluable  in  the  manipulation  and  processing  of  the  experimental  data.  Thanks  to  Mr.  R. 
M.  Harkins  for  his  help  with  LINUS  in  the  lab  and  during  the  deployments  in  the  field. 
Finally,  I  would  like  to  express  my  deep  gratitude  to  my  father  who,  by  his  support  and 
direction  in  my  early  years,  was  one  of  the  primary  reasons  for  my  presence  in  Monterey 
and  subsequently  for  this  thesis. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xviii 


I.  INTRODUCTION 


Advancements  in  semiconductor  technology  during  the  past  two  decades  were  a 
big  impetus  to  Remote  Sensing.  As  sensor  properties  and  com  putational  power  increases. 
Remote  Sensing  is  finding  applications  in  an  increasing  number  of  areas  outside  its 
traditional  boundaries.  The  most  heavily  used  part  of  the  electromagnetic  spectrum  was 
the  visible  region  followed  later  by  the  infrared  and  mic  rowave  regions.  The  ultraviolet 
region,  although  the  least  used  so  far,  gives  the  potential  of  many  new  applications. 

Remote  Sensing  of  the  atmosphere  and  its  constituents  has  long  been  of  great 
importance  for  many  reasons,  both  civilian  and  military.  The  ability  to  detect,  quantify 
and  calculate  the  distribution  of  the  atmospheric  constituents  is  of  great  significance.  In 
the  civilian  case,  monitoring  of  the  atmosphere  is  directly  related  to  public  safety,  quality 
of  life  and  economy,  for  example,  for  weather  prediction,  the  monitoring  of  ozone,  green 
house  effect  gases  and  pollution,  and  the  detection  of  accidents  and  cataclysmic  events 
that  affect  the  atmosphere,  such  as  volcano  eruptions.  In  the  case  of  the  military,  Remote 
Sensing  of  the  atmosphere  can  be  used  for  tactical  and  strategic  intelligence  and 
battlefield  awareness  by  finding  chemical/biological  weapons  gas  clouds  from  safe 
distances,  and  rocket  and  aircraft  exhaust  plumes,  for  example. 

The  Physics  Department  of  the  Naval  Postgraduate  School  has  over  ten  years  of 
experience  with  the  design  and  development  of  numerous  devices  due  to  research 
projects  related  to  atmospheric  Remote  Sensing  in  the  UV.  The  first  device,  beginning  in 
1990,  was  MUSTANG  (Middle  Ultraviolet  SpecTrograph  for  Analysis  of  Nitrogen 
Gases)  and  continued  with  DUUVIS  (Dual  Use  Ultraviolet  Imaging  Spectrometer), 
NUVIS  (Naval  Postgraduate  School  Ultraviolet  Imaging  Spectrometer),  with  the  latest 
being  LINUS  (Lineate  Imaging  Ultraviolet  Spectrometer). 

LINUS  is  an  imaging  spectrometer  with  the  ability  to  operate  from  the  UV 
through  Near  UV  to  the  first  half  of  the  visible  spectrum  depending  on  its  configuration 
and  calibration.  Currently,  it  is  set-up  to  operate  at  approximately  300  ±10  nm  with  the 
purpose  of  detecting  Sulfur  Dioxide  (SO2).  SO2  is  a  trace  atmospheric  gas,  which  can  be 
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found  in  urban  areas  as  a  primary  pollutant  in  the  exhaust  plums  of  industrial 
smokestacks  and  in  volcanic  plumes.  It  has  strong  absorption  bands  centered  at  300nm. 

The  purposes  of  this  thesis  effort  were  to  develop  the  operational  capability  of 
LINUS,  deploy  it  in  the  field  in  order  to  acquire  data  of  natural  plumes,  and  finally 
analyze  the  data  to  specify  the  plume  content  in  SO2.  The  device  was  deployed  at  Lassen 
Volcanic  National  Park  for  experimental  data  acquisition  in  August  and  September  of 
2002.  The  measurements  were  also  modeled  with  the  MODTRAN4  code,  so  that  the 
experimental  data  can  be  better  evaluated  through  analysis  and  direct  comparison. 

This  thesis  consists  of  nine  chapters  and  3  appendices.  Chapter  II  presents  some 
definitions  and  the  basic  physics  involved  in  the  operation  of  LINUS  and  analysis  of  its 
measurements.  Chapter  III  is  dedicated  to  S02,  which  is  the  subject  of  the  measurements. 
It  examines  the  absorption  behavior  of  the  gas  h  the  spectral  region  of  interest  as 
documented  by  detailed  measurements  of  its  absorption  cross-section.  Chapter  IV 
presents  the  current  capabilities  and  practices  worldwide  in  the  detection  and 
measurement  of  volcanic  SO 2  in  the  UV.  Chapter  V  gives  a  brief  description  of  LINUS 
and  its  operation.  Detailed  information  can  be  found  in  other  theses:  Kompatzki,  2000; 
Blackwell,  2001,  Kuriger,  2001  and  Gray,  2002.  Chapter  VI  describes  the  deployment 
and  experimental  data  collection  with  LINUS  at  Lassen  Volcanic  Park  and  also  provides 
a  first  indication  of  the  raw  data.  Chapter  VII  presents  the  concepts  and  technical  details 
of  modeling  the  LINUS  data  with  the  MODTRAN4  code.  The  analysis  of  both  the 
simulated  and  experimental  data  is  presented  in  chapter  VIII.  Finally,  chapter  IX 
summarizes  the  conclusions  deduced  by  this  project. 
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II.  DEFINITIONS  AND  BASIC  PHYSICS 


A.  REMOTE  SENSING 

Remote  Sensing  is  the  acquisition  of  information  about  an  object  without  physical 
contact  (Manual  of  Remote  Sensing,  1983).  The  information  is  acquired  by  detecting  and 
measuring  changes  that  the  object  imposes  on  the  surrounding  field.  This  can  include 
electromagnetic  or  acoustic  energy  emitted  or  reflected  by  the  object  or  changes  to  the 
gravitational  or  magnetic  field  due  to  the  presence  of  the  object  (Elachi,  1987).  Remote 
sensing  measurements  that  depend  on  a  natural  energy  source  are  called  passive,  while 
those  that  depend  on  a  man-made  energy  source,  which  usually  is  on  the  same  platform 
as  the  sensor,  are  called  active. 

Remote  sensing  in  the  electromagnetic  spectrum  depends  upon  observed  spectral 
differences  in  the  energy  reflected  or  emitted  from  features  of  interest.  The  observation  of 
features  at  several  different  wavelengths  in  an  effort  to  derive  information  about  the 
features  and  their  distributions  is  the  basis  of  multispectral  remote  sensing  (Campbell, 
1987). 

With  the  developments  in  detector  technologies,  new  detectors  termed  imaging 
spectrometers  are  available.  They  are  capable  of  recording  image  data  in  a  large  number, 
typically  hundreds,  of  spectral  channels.  For  a  given  pixel,  enough  samples  of  its 
reflectance  properties  can  be  obtained  to  allow  a  very  accurate  characterization  of  the 
pixel’s  spectral  reflectance  or  emittance  curve,  over  the  spectral  region  covered  by  the 
sensor.  This  is  a  case  of  hyperspectral  remote  sensing.  Data  analysis  for  the  hyperspectral 
data  sets  can  extend  beyond  image  analysis  and  involve  procedures  for  the  recognition  of 
fine  spectral  features  in  recorded  spectra  of  individual  pixels.  These  procedures  range 
from  using  models  based  on  a  scientific  understanding  of  the  spectroscopic  processes 
involved  in  spectrum  formation,  to  library  searching  techniques  that  allow  identification 
based  on  matches  to  library  prototypes  (Richards  &  Jia,  1999).  Hyperspectral  imagers, 
such  as  LINUS,  produce  a  spectral  image  of  the  scene  where  each  pixel  has  three 
components:  two  spatial  (x  and  y  coordinates)  and  one  spectral.  Usually  the  sensor 
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records  one  of  the  spatial  components  and  the  spectral  component  and  scans  the  other 
spatial  dimension.  The  data  can  be  thought  as  a  hyperspectral  data  cube. 

Even  higher  spectral  resolution  can  be  achieved  by  ultraspectral  imagers,  which 
can  record  subtle  spectral  details  of  materials.  Figure  2.1  illustrates  the  classification  of 
the  three  remote  sensing  areas. 
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Spectral  discrimination  by  species, 
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Identification  and  discrimination  of  subtle 
spectral  details  of  materials,  vapors  and  aerosols 


Figure  2.1.  Spectral  imager  classifications.  (From:  Multispectral  Users  Guide,  1995) 


B.  ELECTROMAGNETIC  RADIATION 
1.  Electromagnetic  Spectrum 

Although  the  electromagnetic  spectrum  is  continuous  in  nature,  major  divisions 
and  subdivisions  have  been  established  for  convenience  and  by  tradition.  Figure  2.1 
shows  the  typical  division  of  the  spectrum  in  terms  of  frequency,  wavelength  and  photon 
energy.  The  spectral  regions  most  heavily  used  in  remote  sensing  are  the  visible  and 
infrared  but  the  ultraviolet  and  microwave  regions  are  also  used. 
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Figure  2.2.  Electromagnetic  spectrum.  (From:  Barrett  &  Curtis,  1982) 


2.  Blackbody  Radiation 

All  objects  at  temperatures  above  absolute  zero  radiate  energy.  When  they  are  in 
thermal  equilibrium  with  their  surroundings,  they  absorb  and  radiate  the  same  amount  of 
energy  and  their  temperature  stays  constant  with  time.  As  the  temperature  of  the  object 
increases,  the  spectral  distribution  and  amplitude  of  the  radiated  energy  changes.  In  order 
to  approximate  this  behavior,  the  concept  of  blackbody  was  introduced.  Blackbody  is  an 
ideal  source  of  electromagnetic  energy,  with  radiation  properties  that  depend  only  on  its 
temperature.  The  spectral  distribution  of  the  radiated  intensity  R,  is  given  by  Planck’s 
radiation  law: 


R(X)  = 


2k  he2  1 


(2.1) 


where  h  is  the  Planck’s  constant,  c  is  the  speed  of  light,  X  is  the  wavelength,  k  is  the 
Boltzman’s  constant  and  T  is  the  absolute  temperature  of  the  blackbody.  Figure  2.3 
illustrates  the  behavior  of  a  blackbody  at  various  temperatures. 
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Figure  2.3.  Blackbody  radiation  curves  for  various  temperatures.  As  shown  by  (a)  and 
(b)  as  the  temperature  increases  there  is  a  significant  change  of  scale  and  the  maximum 
irradiance  occurs  at  shorter  wavelengths.  (From:  Elachi.  1987) 


The  total  radiated  energy  over  the  entire  spectrum  is  given  by  Stefan  -Boltzmann’s  law: 

R  =  oT4  (2.2) 

where  c  =5.67xl0~8  [W/(m2K4)]  and  the  maximum  of  the  spectral  distribution  is  given  by 
Wien’s  displacement  law: 
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^max^  =  2898  [p,mK] 


(2.3) 


The  concept  of  blackbody  is  important  because  it  c  an  be  used  to  model  the  sun, 
which  is  the  source  for  passive  remote  sensing.  Figure  2.6  shows  that  the  sun  radiates 
approximately  as  a  blackbody  at  about  5900K. 


3.  Interactions  with  the  Atmosphere 

Since  all  radiation  used  for  remote  sensing  of  terrestrial  objects  must  pass  through 
the  atmosphere,  interactions  with  it  are  very  important.  The  atmosphere  consists  of  a 
generally  homogeneous  mixture  of  0>,  N2  and  Ar  with  variable  concentrations  of  water 
vapor,  C02  and  O3.  The  pressure  and  density  of  the  atmosph  ere  decrease  approximately 
exponentially  with  altitude  while  the  temperature  has  a  more  complicated  variation, 
which  characterizes  the  division  of  the  atmosphere  into  layers.  The  composition  of  the 
atmosphere  is  shown  in  table  2.1  and  its  layers  are  shown  in  figure  2.4. 


Constituent  Name  /  (type) 

Formula 

Molec.  Mass 

%  by  volume 

Nitrogen/  (Permanent) 

n2 

28.0134 

78.084 

Oxygen  /  (Permanent) 

o2 

31.9988 

20.9476 

Argon  /  (Permanent) 

Ar 

39.948 

0.934 

Water  Vapor  /  (Variable) 

H20 

18.0160 

0-7 

Carbon  Dioxide  /  (Variable) 

CO  2 

44.00995 

0.01-0.1 

Ozone  /  (Variable) 

o3 

47.9982 

0-0.01 

Neon  /  (Permanent) 

Ne 

20.183 

0.001818 

Helium  /  (Permanent) 

He 

4.0026 

0.000524 

Methane  /  (Permanent) 

ch4 

16.04303 

0.0002 

Sulfur  Dioxide  /  (Variable) 

SO  2 

64.064 

0-0.0001 

Hydrogen  /  (Permanent) 

h2 

2.01594 

0.00005 

Nitrogen  Dioxide  /  (Variable) 

no2 

46.0055 

0-0.000002 

Table  2.1.  Atmospheric  constituents.  (From:  Durkee,  1999) 
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Figure  2.4.  Atmospheric  layers.  (From:  Lenoble,  1993) 

The  main  processes  between  the  electromagnetic  radiation  and  the  atmosphere  are 
graphically  illustrated  in  figure  2.5  below  and  are  as  follows: 


Figure  2.5.  Interaction  of  EM  radiation  with  the  atmosphere.  (From:  Gibson,  2000) 
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a.  Scattering 

Scattering  is  the  redirection  of  the  radiation  by  molecules  of  atmospheric 
gases  or  by  particles  suspended  in  the  atmosphere.  The  amount  of  scattering  depends  on 
the  size  of  the  molecules  or  particles,  their  abundance,  the  wavelength  of  the  radiation 
and  the  depth  of  the  atmosphere  through  which  the  radiation  is  traveling. 

Rayleigh  scattering  occurs  when  molecules  of  atmospheric  gases  or 
particles  have  diameters  that  are  very  small  relative  to  the  wavelength  of  the  radiation.  It 
is  wavelength  dependent  with  a  factor  of  (l/A,4) ,  which  means  that  short  wavelengths  are 

scattered  more.  It  is  the  dominant  scattering  process  above  the  troposphere  and  is 
responsible  for  the  existence  of  the  UV  background  of  the  sky  used  for  LINUS 
measurements. 

Mie  scattering  is  caused  by  larger  particles,  such  as  dust  smoke  and  water 
droplets,  that  their  diameters  are  comparable  to  the  wavelength  of  electromagnetic 
radiation.  It  is  wavelength  dependent  but  in  a  very  complicated  way.  It  is  dominant  in  the 
lower  atmosphere,  where  larger  particles  are  more  abundant. 

Nonselective  scattering  is  caused  by  particles  that  are  much  larger  than  the 
wavelength  of  the  scattered  radiation,  such  as  larger  water  droplets  or  large  particles  of 
dust.  It  is  not  wavelength  dependent  and  causes  a  whitish  or  grayish  haze  at  visible 
wavelengths  (Campbell,  1987). 

b.  Transmission 

Transmission  is  the  process  when  incident  radiation  passes  through  the 
atmosphere  without  measurable  attenuation.  As  the  radiation  passes  through  layers  of  the 
atmosphere  with  different  constituents  and  densities,  it  is  refracted  (deflected).  Refraction 
affects  the  speed  (decreases)  and  the  wavelength  (increases)  of  the  radiation  but  not  its 
frequency. 


c.  Absorption 

Absorption  is  an  important  process  since  it  affects  the  electromagnetic 
radiation  available  for  remote  sensing  and  also  -  for  the  purpose  of  this  thesis  -  provides 
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the  basis  for  detection  and  quantitative  estimation  of  the  S02.  Through  this  process,  the 
radiation  is  absorbed  by  molecules  and  it  raises  the  temperature  of  the  atmosphere,  causes 
chemical  reactions  or  is  reradiated  at  longer  wavelengths.  The  main  molecules 
responsible  for  the  atmospheric  absorption  are  Q,  Q>,  H.0  and  C02.  The  combined 
result  of  scattering  and  absorption  to  the  amount  of  solar  radiation  reaching  the  surface  of 
the  earth  at  each  wavelength  is  illustrated  in  figure  2.6  below. 
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Figure  2.6.  Solar  irradiance  approximated  by  a  blackbody,  at  the  top  of  the 
atmosphere  (TOA)  and  at  sea  level.  (From:  Elachi,  1987) 

d.  Reflection 

Reflection  is  mainly  caused  by  clouds  and  is  not  a  topic  of  inter  est  for  this 

thesis. 
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III.  SULFUR  DIOXIDE 


A.  GENERAL 

Sulfur  dioxide  (SO2)  is  a  trace  chemical  compound  in  the  Earth’s  atmosphere, 
mainly  present  in  the  troposphere  where  it  is  a  primary  pollutant  released  by  fossil  fuel 
combustion,  volcanic  emissions,  oxidation  of  organic  material  in  soils  and  biogenic 
emissions  over  the  oceans.  Some  estimates  consider  that  coal  burning  accounts  for  50% 
of  the  annual  global  emissions  and  oil  burning  for  another  25  -30%.  Volcanic  emissions 
due  to  eruptions  are  responsible  for  the  presence  of  S02  in  the  stratosphere. 

SO2  is  a  colorless  gas  with  a  suffocating,  choking  odor  and  its  density  is  over 
twice  of  that  of  air  at  the  same  conditions  (2.618  g/lt  at  25  °C  and  1  atm).  It  reacts  on  the 
surface  of  various  airborne  solid  particles  (aerosols)  and  can  be  oxidized  within  airborne 
water  droplets,  producing  sulfuric  acid,  which  can  be  transported  by  winds  and  is 
deposited  as  acid  rain. 

The  lifetime  of  S02  in  the  troposphere  is  a  few  days,  while  in  the  stratosphere  it  is 
several  weeks.  This  makes  S02  from  volcanoes  one  of  the  most  important  sources  of 
stratospheric  aerosols.  Its  abundance  in  the  atmosphere  is  highly  variable  above  a  very 
low  background  concentration.  Clean  continental  air  contains  less  than  1  ppb  of  S02. 


B.  ABSORPTI  ON  PROCESSES  IN  THE  VISIBLE  AND  UV 

The  absorption  processes  in  the  visible  and  UV  regions  of  the  spectrum  mainly 
involve  electronic  transitions  of  the  valence  electrons  of  the  atoms.  When  an  electronic 
transition  occurs  in  a  molecule,  the  redistribution  of  the  electrons  changes  the 
electrostatic  forces  operating  on  the  nuclei.  The  molecule  responds  to  the  change  of 
forces  by  bursting  into  vibration.  In  other  words,  an  electronic  transition  is  usually 
accompanied  by  a  vibrational  transition.  The  extra  vibrational  structure  in  a  spectrum  can 
be  resolved  if  the  sample  is  gaseous  but  in  a  liquid  or  solid  sample,  the  lines  merge 
together  and  result  in  a  broad,  almost  featureless  line.  The  vibrational  transitions  that 
accompany  the  electronic  transitions  are  thems  elves  accompanied  by  rotational 
transitions.  As  a  result,  the  electronic  spectra  of  gaseous  samples  can  be  very  complicated 
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(Atkins,  1982).  The  considerable  vibrational  structure  of  gaseous  SO2  at  25°C  is  shown  in 
figure  3.1. 


C.  BEER’S  LAW  -  ABSORPTION  CROSS  SECTION 

When  a  beam  of  monochromatic  radiation  of  initial  intensity  I0  passes  through  a 
material,  absorption  takes  place  and  the  beam  leaving  the  material  has  smaller  intensity,  I. 
The  absorption  of  a  photon  of  the  radiation  corresponds  to  a  change  in  the  internal  state 
of  a  molecule  of  the  material  between  two  energy  levels.  The  energy  of  the  transition  is 
determined  by  the  separation  of  the  energy  levels  between  the  two  states. 

Assume  that  the  material  density  is  c.  Then,  the  absorption  of  radiation  by  a  thin 
layer  of  thickness  dx  at  a  distance  x  from  the  surface  is: 

dl  =  -acldx  (3.1) 

where  o  is  the  absorption  cross  section  of  the  material.  The  absorption  cross  section  is  a 
measure  of  a  molecule's  ability  to  absorb  radiation  at  a  specified  wavelength  and  is 
usually  expressed  in  cm2/particle.  It  depends  on  the  kind  of  molecule,  i.e.,  the  material, 
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and  the  wavelength  of  the  radiation.  If  the  path  length  of  the  radiation  through  the 
material  is  /,  then  by  integrating  equation  (3.1)  the  transmitted  radiation  becomes: 

I  =  I0e^c'  (3.2) 

The  above  equation  is  Beer’s  law  and  applies  to  optically  thin  samples.  Another  version 
of  the  law  contains  the  molar  absorption  coefficient  e . 

Beer’s  law  is  the  main  theoretical  tool  used  to  calculate  absorption  cross  sections 
for  various  materials.  In  the  experimental  setups,  a  monochromatic  light  of  known 
intensity  is  directed  into  a  sample  with  specific  characteristics,  such  as  temperature  and 
pressure,  the  outcoming  light  intensity  is  measured,  and  thus,  the  absorption  cross  section 
can  be  calculated. 

D.  S02  ABSORPTION  CROSS  SECTIONS  IN  THE  NEAR  UV 

The  absorption  spectrum  of  S02  between  285  and  315nm  displays  a  series  of 
well-defined  peaks  and  troughs  spaced  approximately  2nm  apart.  This  constitutes  a 
strong  spectral  signature  useful  in  remote  sensing  applications.  For  this  reason,  LINUS 
was  designed  to  operate  in  the  region  290-310nm  with  a  maximum  spectral  response  at 
about  300nm. 

One  application  is  the  Differential  Absorption  Lidar  (DIAL),  which  uses  the 
difference  in  absorption  cross  section  between  an  absorption  peak  and  an  adjacent  trough 
to  measure  S02  concentrations  remotely  for  environmental  monitoring  purposes.  The 
maximum  differential  absorption  cross  section  occurs  close  to  300  nm.  This,  considering 
that  most  DIAL  systems  use  dye  lasers  that  work  most  efficiently  at  this  wavelength, 
raised  the  need  for  very  accurate  measurements  of  the  S02  absorption  cross  section  in  this 
wavelength  region  (Brassington,  1981).  Another  application  is  the  processing  of  the  data 
produced  by  the  TOMS  instrument,  which  will  be  discussed  later  in  this  thesis. 

There  were  several  measurements  reviewed  for  this  region  of  wavelengths  with 
spectral  resolutions  varying  from  moderate  ( ~  0.02nm)  to  high  ( ~  0.002nm).  In  most  of 
these  cases,  the  experimental  setup  included  frequency  doubled  dye  lasers  and  used  cells 
containing  S02  at  atmospheric  temperatures  (~  20°C)  and  pressures.  Measurements  at 
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high  (~  100°C)  and  low  (~6()"C)  temperatures  and  low  pressures  were  also  conducted.  In 
some  cases,  dry  nitrogen  was  used  as  a  buffer  gas  to  simulate  atmospheric  pressure 
broadening,  which  smears  the  fine  spectral  structure  of  pure  SCU- 

The  results  of  some  of  these  measurements  are  presented  in  figures  3.2  through 
3.7  below.  Figure  3.2,  from  the  work  of  Brassington  in  1981,  shows  the  absorption  cross 
section  of  S02  in  the  spectral  region  from  290  to  317nm,  while  figures  3.3  and  3.4 
compare  the  results  of  that  work  with  previous  results  for  the  peaks  at  296,  298  and 
300nm.  The  presented  Brassington  measurements  were  conducted  over  a  broad  spectral 
region  compared  to  previous  works  at  moderate  resolution  (0.01  nm),  low  pressures  (2.3- 
17.6mbars)  and  ambient  temperatures  (292-300K),  without  the  use  of  buffer  gas.  In  the 
work  of  Woods  et  al.,  1980,  high -resolution  (0.002nm)  measurements  were  conducted  in 
the  spectral  region  from  297.0  to  301.0nm  at  20°C  with  the  use  of  nitrogen  buffer  gas  at 
latm  in  order  to  simulate  in  the  laboratory,  conditions  encountered  in  the  lower 
atmosphere  (pressure  broadening).  Figure  3.5  from  that  work  shows  the  absorption  cross 
section  peaks  at  298  and  300nm  in  detail,  while  figure  3.6  shows  the  same  two  peaks  with 
the  in-between  trough  and  with  nitrogen  buffer  gas.  Finally,  figure  3.7  from  Thomas  et 
al.,  shows  the  S02  absorption  cross  section  in  the  region  from  300.0  to  320.4nm  for  room 
and  low  temperature  (295  and  21  OK)  at  moderate  resolution  (0.03nm). 

LINUS,  used  for  the  measurements  in  this  thesis,  is  currently  set  up  to  operate  in 
the  spectral  region  from  290  to  310nm  with  a  maximum  spectral  response  at  about 
300nm.  The  above-mentioned  measurements,  illustrated  in  the  following  figures,  are  of 
great  significance  since  they  provide  the  needed  spectral  signatures  that  will  allow 
recognition  of  potential  S02  detection  by  LINUS  and  subsequent  quantitative 
estimations. 
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Figure  3.3.  Comparison  of  measurements  between  Brassington  (solid  line)  and 
Thomson  et  al.  (dotted  line).  (From:  Brassington,  1981) 
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Figure  3.4.  Comparison  of  measurements  between  Brassington,  Thomson  et  al.  and 
Woods  et  al.  for  the  300.03nm  absorption  peak.  (From:  Brassington,  1981) 
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Figure  3.5.  High  precision  SO2  absorption  cross-section  peaks  at  298  and  300nm  at 
20°C.  The  dotted  line  is  with  l.OlxltF  Pa  of  nitrogen  buffer  gas.  The  solid  lines  are  is  for 
SO2  alone.  The  original  article  refers  to  the  solid  line  as  the  one  with  the  buffer  gas  and 
the  dotted  as  SO2  alone,  but  this  is  obviously  opposite  since  the  buffer  gas  smoothes  out 
the  spectral  details.  (From:  Woods  et  al.,  1980). 
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Figure  3.6.  High  resolution  absorption  cross  section  of  SO2  at  20°C  with  nitrogen 
buffer  gas  over  the  spectral  region  297.8-300.3nm.  (From:  Woods  et  al.,  1980) 
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IV.  UV  REMOTE  SENSING  OF  VOLCANIC  S02 


A.  GENERAL 

Volcanic  gases  are  a  fundamental  pail  of  volcanic  activity.  Ranging  from 
dramatic  emissions  of  acidic  vapors  released  by  highly  active  volcanoes  to  invisible 
emanations  through  soils,  volcanic  gases  play  an  important  role  in  influencing  the 
behavior  of  a  volcano.  Quantitative  information  on  volcanic  emissions  is  essential  to 
studies  of  climate  change,  volcanic  processes  and  natural  hazards. 

Each  year  volcanic  eruptions  inject  significant  but  variable  amounts  of  ash  and 
gases  into  the  atmosphere,  which  may  encircle  the  globe  for  weeks  after  an  eruption 
endangering  aircraft  and  affecting  the  weather.  The  conversion  of  SCE  to  sulfate  aerosols 
impacts  stratospheric  chemistry  and  the  global  climate.  Volcanoes  are  by  far  the  lar  gest 
contributor  of  stratospheric  sulfate  aerosols.  The  presence  of  aerosols  increases  the 
atmospheric  albedo,  reflecting  solar  radiation,  absorbing  terrestrial  infrared  radiation  and 
modifying  the  greenhouse  effect  in  the  atmosphere.  Aviation  is  also  affected  by  these 
emissions  because  jet  engines  and  avionic  systems  fail  in  flight  when  volcanic  ash  is 
ingested  and  longer  term  damage  results  from  exposure  to  sulfuric  acid  aerosols, 
including  the  crazing  of  aircraft  windows,  fading  of  polyurethane  paint  and  accumulation 
of  sulfate  deposits  in  jet  engines.  By  studying  gases,  volcanologists  gain  insight  into  a 
number  of  volcanic  processes,  and  by  studying  long  -term  gas  data  along  with  other  types 
of  data,  try  to  forecast  eruptive  activity. 


B.  TOMS 

The  presence  of  ozone  (O3)  in  the  earth’s  atmosphere  was  originally  postulated  to 
explain  the  cutoff  of  solar  radiation  reaching  the  ground  at  UV  wavelengths. 
Subsequently  the  region  from  300  to  340nm,  where  strong  ozone  spectral  signature 
exists,  became  the  standard  for  measurements  of  column  ozone  amounts  with  ground  and 
spaceborne  instruments  (Krueger  et  al.,  2000).  The  best  long-term  global  dataset  of 
atmospheric  ozone  was  collected  by  the  Total  Ozone  Mapping  Spectrometer  (TOMS) 
onboard  the  Nimbus -7  satellite.  TOMS  provided  global  ozone  data  almost  every  day 


19 


from  October  1978  to  May  1993  from  a  low  earth  polar  orbit.  The  instrument  operated  by 
comparing  incident  solar  radiation  with  radiation  reflected  from  the  earth’s  atmosphere. 
Six  UV  wavelengths  were  observed,  each  with  a  bandwidth  of  lnm:  313,  318,  331,  340, 
360  and  380nm  (http://daac.gsfc.nasa.gov). 

UV  remote  sensing  of  volcanic  clouds  from  space  began  with  the  discovery  that 
SO2,  in  the  1982  El  Chichon  eruption  clouds,  could  be  mapped  with  TOMS.  During  the 
eruption,  ozone  levels  were  observed  to  be  anomalously  high,  until  it  was  realized  by 
Arlin  Krueger  and  his  colleagues  at  NASA,  that  TOMS  was  detecting  volcanic  SO2  in 
addition  to  ozone.  Ozone  and  SO2  have  similar  absorption  spectra  in  parts  of  the 
wavelength  region  where  TOMS  operates,  as  shown  in  figure  4.1,  with  SO2  producing 
stronger  absorption  at  the  shortest  TOMS  wavelengths  (300-320nm). 


PHCrOABSGRnriON  CROSS  SECTION- 


Figure  4.1.  Absorption  cross  sections  of  atmospheric  UV  absorbers.  (From: 

Strickland,  1999) 


Normally  the  background  amounts  of  SO2  are  so  small  that  they  can  be  neglected 
in  ozone  data  processing.  Flowever,  the  great  abundance  of  SO 2  in  explosive  eruption 
clouds  can  overwhelm  the  absorption  by  ozone  causing  large  errors  in  the  ozone 
retrievals.  This  was  first  observed  in  1982  as  mentioned  previously,  and  led  to  the 
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development  of  new  algorithms  to  separately  determine  ozone  and  S02  in  volcanic 
clouds.  The  spatial  resolution  of  TOMS  is  50km  at  nadir.  Thus,  volcanic  clouds  that 
extend  from  less  than  50km  to  many  thousands  of  kilometers  can  be  measured  (Krueger 
et  al.,  2000).  An  example  of  S02  detection  by  TOMS  shown  in  figure  4.2  is  the  volcanic 
cloud  created  after  the  eruption  of  Mt.  Pinatubo  in  the  Philippines.  The  units  in  the  figure 
are  Dobson  Units  (DU):  lDU=lmatm-cm=  lOppm-m. 
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Figure  4.2.  False-color  image  of  S02  from  the  June  16,  1991  eruption  of  Mt. 
Pinatubo,  Philippines,  produced  from  TOMS  data.  The  triangular  symbol  marks  the 
location  of  the  source  volcano.  (From:  <http://skye.gsfc  .nasa.gov>,  2002) 


Table  4.1  below  shows  the  satellite  missions  with  TOMS  instruments.  Earth 
Probe  TOMS  is  currently  the  only  NASA  spacecraft  in  orbit  specializing  in  ozone 
retrieval. 
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Satellite 

Launch/ 

Activation  Dates 

Ending  Date 

Orbit  Height,  km 

LT  Ascending 
Node 

Nadir 

Footprint,  km 

Nimbus-7 

October  24.  1978 

May  6.  1993 

955 

1 1  JO  AM 

50 

November  1.  1978 

Mcttor-3/5 

Aug  IS.  1991 

December  24,  1994  1200 

Not  Sun- 

62 

Aug  22. 1991 

synchronous 

TOMS  Earth 

July  2.  1996 

Operating, 

500 (Launch- 

11:30  AM 

24 

Probe 

July  17.  1996 

July  1999 

September  1997) 

739  (November  1997 

-present) 

39 

ADEOS 

August  17.  1996 

June  30.  1997 

797 

10:30  EM 

42 

September  11. 1996 

TBD  platform 

2002 

Table  4.1.  Satellite  missions  with  TOMS  instruments.  (From:  Krueger  et  al.,  2000) 

C.  BREWER  SPECTROMETER 

The  Brewer  spectrometer  (see  figure  4.3)  is  a  ground  based  UV  spectrometer 
designed  for  automated  monitoring  of  total  ozone  and  SO 2  by  measuring  the  spectral 
absorption  of  direct  sunlight  in  five  narrow  bands  between  305  and  320nm.  In  addition, 
the  instrument  measures  direct  and  global  (direct  +  diffuse)  spectral  irradiance  in  the 
wavelength  range  from  290  to  372nm.  If  direct  sun  is  available,  total  ozone  is  derived 
from  direct  sun  measurements.  For  cloudy  sky  conditio  ns,  ozone  may  be  derived  from 
scattered  UV  radiation  from  the  zenith  direction.  The  center  wavelengths  for  the  5  UV 
channels  are  306.3,  310.1,  313.5,  316.8  and  320.1nm. 

Successful  measurement  of  column  SO  2  amounts  in  volcanic  clouds  depends  on  a 
chance  passage  of  the  cloud  over  a  station.  One  of  the  few  cases  in  which  this  occurred 
was  when  the  Mt.  Spun-,  Alaska  cloud  crossed  over  Toronto  on  September  19,  1992.  A 
Brewer  spectrometer  was  looking  up  at  the  plume  and  allowed  both  the  Nimbus  -7  and 
Meteor- 3  TOMS  instruments  to  be  directly  compared  with  ground  observations. 

An  example  of  the  output  produced  by  the  Brewer  station  of  Montreal,  Canada  is 
shown  in  figure  4.4.  The  graph  shows  daily  mean  total  ozone  values  (in  Dobson  Units)  as 
a  function  of  the  season.  The  results  of  the  most  accurate  Brewer  direct  sun 
measurements  are  shown  in  red.  All  other  types  of  data  are  shown  in  blue.  The  green  line 
indicates  total  ozone  values  smoothed  by  a  one  month  running  mean.  The  dashed  lines 
represent  the  “normal”,  pre-1980  total  ozone  annual  cycle  and  two  standard  deviations 
envelope  for  daily  mean  total  ozone  (From  the  World  Ozone  and  Ultraviolet  Radiation 
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Data  Center  (WOUDC)).  If  a  volcanic  SO 2  cloud  happens  to  pass  over  the  station,  these 
data  are  modified  due  to  the  additional  absorption  and  S02  quantitative  calculations  can 
be  made. 


Figure  4.3.  Brewer  Station  on  the  roof  of  the  Environment  Canada  Building,  Toronto. 
(From:  <http://www.ns.ec.gc.ca>,  2002) 
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Figure  4.4.  Example  of  Brewer  station  output.  (From:  WOUDC,  2002) 
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D.  COSPEC 


COSPEC  (see  figure  4.5)  is  an  instrument  developed  in  the  1960s  by  Barringer 
Research,  a  company  in  Toronto,  Canada,  to  measure  emissions  of  SO  2  from  industrial 
stacks  for  pollution  monitoring  purposes.  Volcanologists  quickly  realized  the  potential  of 
COSPEC  since  volcanoes  are  essentially  natural  smokestacks.  It  was  first  used  in  1971  at 
Mt.  Mihara  in  Japan  and  then  usage  spread  around  the  world  providing  a  collection  of 
data  that  has  been  used  to  estimate  a  worldwide  SO2  flux  (emission  rate)  from  volcanoes. 


Figure  4.5.  Close-up  view  of  COSPEC.  (From:  USGS  website,  2002) 

1.  Brief  Description  of  Instrument  and  Its  Operation 

The  instrument  operates  in  the  UV  region  of  the  spectrum  using  scattered  UV 
radiation  from  the  sky.  The  UV  sensor  has  an  effective  spectral  range  of  280-320nm  but 
the  instrument  is  usually  operated  from  300  to  315nm.  Its  operation  is  based  on  the 
principles  of  correlation  spectroscopy  and  its  internal  design  is  illustrated  in  figure  4.6. 
Radiation  is  gathered  by  the  telescope  with  a  field  of  view  of  23  milliradians  by  7 
milliradians  (1.318  by  0.401  degrees).  The  collected  radiation  is  focused  on  a  diffraction 
grating,  which  separates  the  radiation  into  individual  wavelengths.  This  radiation  is 
focused  on  the  rotating  correlator  disk,  which  has  four  sets  of  seven  narrow  slits  (four 
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masks)  engraved  on  it.  The  rotation  sequentially  places  the  slits  so  that  they  coincide 
approximately  with  the  peaks  and  troughs  of  the  S02  absorption  bands.  This  modulated 
radiation  is  monitored  by  the  photomultiplier  tube  and  the  resulting  voltages  are 
processed  by  the  electronics.  The  ratio  between  the  two  sets  of  incoming  radiation  (S02 
peaks  and  troughs)  is  known  when  no  S02  is  present  to  absorb  energy.  When  SO  2  is 


Figure  4.6.  Internal  design  of  COSPEC.  (From:  COSPEC  IV  Manual,  1976) 

present,  one  set  of  the  slits  (corresponding  to  the  troughs)  loses  ener  gy,  and  thus,  the  ratio 
changes.  These  changes  are  electronically  processed  to  yield  an  output  voltage,  which  is  a 
function  of  the  concentration  times  the  pathlength  (optical  depth)  of  the  S02  within  the 
field  of  view  of  the  device  (Millan  et  al.,  1976  and  Stoiber  et  al.,  1983).  The  instrument  is 
connected  to  a  data  logger  or  portable  computer,  which  records  the  output  in  digital 
format  and  an  analogue  paper  chart  recorder,  which  plots  the  voltage  recordings  that 
correspond  to  a  real  time  S02  profile  for  the  area  covered  by  the  sensor.  An  example  of 
the  paper  chart  recorder  output  is  shown  in  figure  4.7.  The  instrument  response  is 
calibrated  by  the  rotation  of  cells,  which  contain  known  S02  concentrations,  into  the 
optical  path.  Frequent  calibrations  indicate  the  amount  of  instrument  drift. 
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Figure  4.7  Example  of  COSPEC  output  to  a  paper  chart  recorder.  (From: 
<http://www.gsj.jp>,  2002) 


2.  Calculations  and  Limitations  of  the  SO  2  Flux  Determination 

The  initial  S02  absorption  profile  is  converted  from  volts  to  an  average 
concentration  path  length  in  parts  per  million-meter  (ppm-m).  One  ppm-m  of  S02  is 
equivalent  to  one  cubic  centimeter  of  S02  gas  uniformly  mixed  in  one  million  cubic 
centimeters  of  air  that  is  viewed  by  COSPEC  over  an  optical  path  of  one  meter  at  a 
pressure  of  101.325  kPa  and  a  temperature  of  20°C  (Jones  and  Stix,  2001).  The 
concentration  is  multiplied  by  the  width  of  the  plume  cross  -section  and  produces  a  S02 
cross-section  in  ppm-m2.  The  cross-section  is  then  multiplied  by  the  speed  at  which  the 
plume  crosses  this  two  dimensional  area  (usually  the  wind  speed  is  used)  producing  the 
S02  flux  in  ppm-m3/sec.  Finally  this  value  is  converted  to  metric  tons  per  day  (t/d)  for 
convenience  in  reporting  (Stoiber  et  al.,  1983).  Detailed  information  about  the 
calculations  can  be  found  among  other  sources  in  Mares,  [2002].  The  S02  flux 
calculation  is  summarized  in  table  4.2  and  an  example  of  the  final  product,  -  S02 
emission  rate  in  tons  per  day  for  the  Kilauea  Summit,  Hawaii  -  is  shown  in  figure  4.8 
below. 
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COSPKC  Output  :  ppm. n 
Pluxie  Width  :  m 

Plurx*  Speed  :  m .  as 

COSPEC  Output  x  Width  of  Pluxe  “  SO-  Plume  Cross-Section 

2 

ppro.m  x  m  -  ppm. n 

.  Pluttti  Cross-Section  x  Plune  Spued  •»  SO-  Flux 

2-13-1 

pprr.n  x  i».s  -  ppm.m  .  * 

Conversion  to  metric  con*. day  '  : 

(D«nsity  of  SO^  gas  :  2.86  x  10  ^  g.ppm  Ex'1  ;  STP  correction  factor  :  273/293) 
X  ppm. ra^  s  1  x  2.86x10  3  g.ppm  1  *  x  10  6  tona.R  1  x  86,400  s.duy*'  x  273.293 

•  Y  tons. day  * 


Table  4.2.  Calculation  of  S02  flux  from  COSPEC  output.  (From:  Stoiber  et  al.,  1983) 
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Figure  4.8.  Example  of  long-term  COSPEC  calculation  results  in  metric  tons/day 
(mt/d)  of  S02.  Kilauea  Summit,  Hawaii  S02  emission  rates.  (From:  Elias  et  al.,  1998) 

The  estimation  of  the  error  in  the  S02  flux  calculations  by  Stoiber  et  al.  is 
presented  in  table  4.3.  An  example  of  flux  calculation  together  with  estimated 
uncertainties  for  Kilauea’ s  east  rift  zone,  Hawaii  is  shown  in  figure  4.9. 
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Daca  : 


Calibration  Cell  Concantrations  :  425  t  12  ppn.m 
Chart  Record  Reading  Error  i  t  0.5  nn 

(0.5  mm  *  6  ppm. a) 

For  an  average  deflection  of  IOO  ppn.ei 
Distance  Determination  Error 
Wind speed  determination 
(generally) 

(worst  case) 


13  1 


1  6  X 
t  5  *  to  1  10  X 

110  t  to  i  20  1 
1  40  X 


Total  Error  >  (square  root  of  Che  sum  of  the  squares) 

Generally  :  •  r13  X  to  t  23  X 

Worst  Case  :  ■  *  42  X 


Table  4.3.  Limitations  of  the  COSPEC  calculations.  (From:  Stoiber  et  al.,  1983) 


Figure  4.9.  Example  of  estimated  uncertainties  of  SCL  flux  calculations  based  on 
vehicle -mounted  COSPEC.  The  black  vertical  bars  represent  the  standard  deviation  of  all 
traverses  on  a  single  day.  Kilauea’s  East  Rift  Zone,  Hawaii.  (From:  Elias  et  al.,  1998) 

The  uncertainty  of  wind  speed  is  the  largest  source  of  error,  which  increases  with 
the  variability  of  the  wind  at  low  wind  speeds.  Other  factors  affecting  the  quality  of  the 
COSPEC  data  and  calculations  (Millan  et  al.,  1976  and  Stoiber  et  al.,  1983)  are  the 
following: 

•  Due  to  weather  conditions  (mainly  shear  winds)  the  plume  profile  may  be 
different  than  the  Lagrangian  and  Eulerian  profiles  assumed  in  order  to 
derive  the  methodology  used  for  the  calculations. 
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•  The  time  of  day  that  the  measurements  are  taken.  Since  scattered  solar 
radiation  is  used  as  the  source,  early  or  late  in  the  day  measurements  are 
affected  by  the  lower  levels  of  UV  radiation,  which  result  in  a  decreased 
signal  to  noise  ratio,  and  thus,  the  precision  of  the  measurements  is 
reduced. 

•  The  lack  of  light  due  to  weather  conditions  (clouds,  fog,  etc.). 

•  The  presence  of  volcanic  ash  in  the  plume  does  not  allow  complete 
transmission  of  UV  radiation,  resulting  in  a  reduction  in  the  measur  ed  SO2 
flux  from  the  actual  value. 

•  The  target  gas  signal  is  diluted  by  the  radiation  scattered  by  the 
atmosphere  and  condensed  water  vapor  into  the  field  of  view  from  outside 
the  solid  angle,  which  collects  radiation  observed  by  the  sensor.  This 
effect  depends  on  the  distance  of  the  instrument  from  the  plume  and  the 
width  of  the  plume.  For  example,  when  one  views  a  vertical  plume,  as 
from  a  stack  or  volcano  with  the  COSPEC  at  a  ground  point  some  distance 
away,  the  signal  may  decrease  as  much  as  60%  at  a  dis  tance  of  only  500m 
from  the  crater,  if  the  SO2  concentration  path-length  is  lOOOppm-m  and 
the  meteorological  range  is  9.8km  (Stoiber  et  al.,  1983).  If  a  plume  is 
viewed  from  below  and  it  is  wider  than  the  vertical  distance  between  the 
sensor  and  the  bottom  of  the  plume,  all  the  UV  radiation  reaching  the 
COSPEC  comes  through  the  plume.  Except  for  the  edges  of  the  plume, 
this  radiation  is  not  scattered  into  the  field  of  view  of  the  sensor  from  the 
outside. 

3.  S02  Flux  Measurement  Techniques  Using  COSPEC 

There  are  three  standard  techniques  of  collecting  data  on  volcanic  plumes  with 
COSPEC:  airborne,  ground  mobile,  and  ground  stationary  (see  figures  4.10  and  4.11). 
COSPEC  has  also  been  used  on  a  boat.  Each  technique  has  advantages  and  disadvantages 
depending  on  the  location,  access  to  it,  surrounding  terrain,  plume  size,  shape  and 
altitude,  cost  and  availability  of  airplanes  etc.  In  the  first  two  techniques,  COSPEC  is 
mounted  on  an  airplane  or  vehicle  that  makes  traverses  beneath  the  gas  plum,  with  the 
traverses  oriented  perpendicular  to  the  direction  of  the  plume.  The  preferred  technique  is 
with  an  airplane  equipped  with  a  good  navigation  system,  since  much  data  can  be 
obtained  in  a  short  period  of  time  and  plumes  that  cannot  be  accessed  by  other  means  can 
be  measured  (Stoiber  et  al.,  1983). 
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Figure  4.10.  COSPEC  mounted  on  vehicle.  (From:  <http://www.volcanoes.ca>,  2002) 


Figure  4.1 1.  COSPEC  on  airplane  and  on  ground.  (From:  USGS  website,  2002) 


The  ground  stationary  technique  will  be  presented  here,  since  it  is  similar  in 
concept  with  the  one  used  with  FINUS.  It  involves  mounting  COSPEC  on  a  tripod  either 
directly  under  or  off  to  the  side  of  the  plume.  If  COSPEC  is  under  the  plume,  then  the 
instrument  is  aligned  with  the  plume  axis  and  a  right -angle  mirror/rotator  is  attached  to 
the  telescope.  The  mirror  is  then  rotated  through  a  series  of  angles  until  the  whole  width 
of  the  plume  has  been  scanned.  When  the  COSPEC  is  to  the  side  of  the  plume,  the 
goniometer  head  on  the  tripod  can  be  used  to  scan  the  COSPEC  through  a  vertical  angle 
if  the  plume  is  horizontal  or  a  horizontal  angle  if  the  plume  is  vertical  as  it  leaves  the 
volcano  (Stoiber  et  al.,  1983). 
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The  ground  stationary  technique  probably  produces  the  least  reliable  data  because 
varying  plume  geometries  often  make  calculations  difficult  and  large  and  variable  path- 
lengths  from  the  plume  cause  variable  signal  attenuation  due  to  scattering,  thus  affecting 
the  measured  SO2  flux.  However,  this  technique  is  the  only  one,  which  can  be  applied 
under  certain  circumstances,  such  as  at  volcanoes  where  there  is  no  effective  air  or 
ground  access  (Stoiber  et  al.,  1983).  A  comparison  between  the  vehicle  and  tripod -based 
measurements  of  the  Kilauea  Volcano,  Hawaii  from  1995  to  1997  showed  that  the 
vehicle -based  measurements  were  1.3-2  times  greater  than  the  tripod -based.  That 
discrepancy  was  attributed  to  a  combination  of  factors  including  multiple  degasing 
sources,  plume  thickness  and  geometry,  distance  from  the  measurement  site  to  the  plume 
and  instrument  configuration.  It  was  concluded  that  the  tripod-based  data,  although  they 
provide  important  information  for  the  vent  activity,  should  be  interpreted  with  care  and 
the  best  estimate  of  SO2  emission  rates  was  provided  by  the  vehicle -based  technique 
(Elias  et  al.,  1998). 

E.  FLYSPEC 

FLYSPEC  is  a  sensor  currently  under  development  by  Keith  Horton,  Assistant 
researcher  at  the  Hawaii  Institute  of  Geophysics  and  Planetology.  It  is  essentially  a 
miniature  COSPEC  (its  name  emphasizes  its  size  in  comparison  to  COSPEC)  that  takes 
advantage  of  new  spectrometer  and  computing  technology.  Apart  from  the  reduced  size 
of  FLYSPEC,  other  differences  with  COSPEC  are  the  lower  cost  of  the  device,  reduced 
data  processing  time  and  real  data  analysis  capability.  Although  no  formal  literature  is 
currently  published,  details  about  the  design  and  operation  of  the  device  can  be  found  in 
Mares,  2002.  Figure  4.12  shows  COSPEC  and  FLYSPEC  mounted  on  the  same  vehicle, 
figure  4.13  shows  some  sample  data  collected  by  FLYSPEC  at  the  Kilauea  volcano, 
Hawaii  in  March  2002  and  figure  4.14  shows  a  comparison  between  COSPEC  and 
FLYSPEC  data  taken  concurrently  at  the  same  site. 
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Figure  4.12.  Concurrent  data  collection  with  COSPEC  and  FLYSPEC.  (From:  Mares, 

2002) 


Figure  4.13.  FFYSPEC  sample  data  frame  collected  at  Kilauea  Volcano,  Hawaii  in 
March  2002.  (From:  Mares,  2002) 
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Figure  4.14.  Comparison  between  COSPEC  and  FLYSPEC  data  taken  concurrently  in 
March  2002  at  Kilauea  Volcano,  Hawaii.  (From:  Mares,  2002) 

F.  SUMMARY 

TOMS  is  a  spacebome  instrument  that  allows  detection  and  quantification  of 
large-scale  SO2  emissions  during  volcanic  eruptions.  The  results  depend  on  the 
algorithms  used  to  extract  the  SO2  information  from  the  raw  data.  TOMS  SO2  data  led  to 
significant  advances  from  determining  eruptive  emission  rates  to  tracking  of  individual 
clouds  and  assessing  global  volcanism  and  atmospheric  impacts. 

The  ground-based  instruments  can  only  locally  detect  and  quantify  S02  emissions 
and  their  accuracy  greatly  depends  on  weather  conditions  and  especially  wind  direction 
and  velocity.  They  use  scattered  UV  radiation  except  for  the  Brewer  spectrometer,  which 
uses  direct  sunlight  when  there  are  no  clouds.  Brewer  spectrometers  are  stationary 
instruments  mainly  used  to  measure  ozone  in  the  atmosphere  and  they  can  only  measure 
volcanic  SO  2  by  chance  if  the  weather  conditions  bring  a  volcanic  cloud  over  the  Brewer 
station.  COSPEC  and  FLYSPEC  are  mobile  instruments,  which  can  be  mounted  on 
vehicles  or  aircraft  that  have  to  pass  under  the  plume,  or  they  can  be  used  stationary  on 
the  ground  at  a  site.  COSPEC  has  a  30-year  history  of  SO2  observations,  has  produced  a 
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large  collection  of  data  and  is  the  main  ground  based  instrument  in  use  today.  FLYSPEC 
is  a  new  instrument  still  under  development,  based  on  the  same  principles  as  COSPEC 
but  takes  advantage  of  new  technologies  in  order  to  be  smaller,  faster  and  cheaper  than  its 
predecessor.  Figure  4.15  below  summarizes  the  data  collection  techniques  for  UV  remote 
sensors  while  figure  4.16  shows  the  wavelength  region  in  which  they  operate. 


OBSERVATIONS  OF  VOLCANIC  S02  CLOUDS 


Figure  4.15.  UV  sensors  used  in  the  detection  of  volcanic  SCE.  (From:  Krueger  at  al., 

2000) 
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Figure  4.16.  Absorption  spectra  of  S02  and  03  and  wavelength  bands  used  by  the  UV 
detectors.  (From:  Krueger  at  al.,  2000) 
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V.  LINEATE  IMAGING  NEAR  ULTRA  VIOLET 
SPECTROMETER  (LINUS) 


A.  GENERAL 

LINUS  is  a  hyperspectral  ultraviolet  imaging  spectrometer  designed  to  study 
atmospheric  gas  plumes.  It  was  designed  and  developed  by  Professor  Scott  Davis  and  his 
students  of  the  Physics  Department  at  the  Naval  Postgraduate  School. 


B.  DEVICE  DESCRIPTION  AND  OPERATION 

In  order  to  understand  the  capabilities  of  LINUS  and  the  kind  of  data  it  can 
produce,  a  brief  description  of  the  device  and  its  operation  is  given  below.  Detailed 
information  about  the  device  can  be  found  in  (Kompatzki,  2000)  and  (Gray,  2002). 
Elements  from  the  first  work  are  shown  below.  The  architecture  of  LINUS  is  shown  in 
figure  5.1. 


1.  Hardware 

LINUS  (see  figure  5.4)  consists  of  four  main  subsystems: 
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a.  Host  Computer 

The  host  computer  is  the  heart  of  LINUS  acting  as  the  main  processing, 
synchronizing  and  control  unit.  It  is  an  ICS  Advent  system  composed  of  a  7490  series 
chassis,  a  14.1”  liquid  crystal  display,  keyboard  and  mouse,  all  ruggetized  for  field  use. 
The  CPU  board  is  based  on  an  Intel  Pentium  processor. 

b.  Scanning  system 

The  scanning  system  performs  the  optomechanical  scanning  of  the  scene. 
It  is  a  closed  loop  angular  position  control  system  that  points  the  image -scanning  mirror 
to  any  desirable  heading.  It  consists  of  the  following: 

•  A  National  Instruments  Flexmotion  6C  motion  card. 

•  A  12FS  model  Kollmorgen  Motion  Technologies  Group  high  torque  low  - 
inertia  DC  servomotor. 

•  A  R137  model  Gurley  Precision  Instruments  digital  encoder,  which 
provides  position  feedback.  The  encoder  has  a  resolution  of  144,000 
counts  per  revolution. 

•  A  National  Instruments  nuDrive  multiaxis  power  amplifier. 

•  The  scanner  mirror  assembly  which  houses  the  mirror,  the  DC  motor  and 
the  encoder. 

c.  Camera  system 

The  camera  is  a  Princeston  Instruments  Intensified  Pentamax  System, 
which  consists  of  an  ultraviolet  intensified  CCD  camera,  a  power  supply,  a  PG-200 
pulser  with  an  MCP-100  high  voltage  power  supply  and  a  high-speed  serial  link  board. 
At  the  front  end  of  the  camera  is  an  integrated,  ultraviolet  enhanced  microchannel  plate 
(MPC)  image  intensifier.  The  output  of  the  intensifier  is  optically  connected  to  a  512x512 
charge  coupled  device  (CCD)  detector  ar  ray  via  a  tapered  fiber  optic  bundle.  The  camera 
is  cooled  to  a  temperature  set  by  the  user  and  its  spectral  range  is  200-550nm.  Camera 
control  and  data  acquisition  is  accomplished  through  the  host  computer,  which  is 
connected  to  the  camera  through  the  hig  h-speed  serial  link. 
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d.  Optics 

The  optical  layout  of  LINUS  is  shown  in  figure  5.2.  At  the  front  end  is  the 
scanner  mirror,  which  is  operated  by  the  servo  system.  The  radiation  reflected  by  the 
mirror  passes  through  a  filter,  which  in  the  case  of  the  measurements  used  in  support  for 
this  thesis,  is  centered  at  about  300nm  with  a  spectral  response  as  shown  in  figure  5.3. 
After  the  filter,  the  primary  telescope  objective  focuses  the  radiation  onto  a  vertical  slit 
aperture.  Then,  a  collimating  lens  sends  the  radiation  to  a  reflective  plane  diffraction 
grating,  and  after  been  focused  by  the  camera  objective,  ends  up  onto  the  input  plane  of 
the  camera’s  image  intensifier. 


Figure  5.2.  LINUS  optical  layout.  (From:  Kompatzki,  2000) 
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2.  Software 

All  the  software  is  installed  on  the  host  computer  and  consists  of  the  following 
main  modules: 

•  Windows  NT  4.0  is  the  operating  system. 

•  Lab  VIEW  5.1  is  used  to  control  the  scanning  subsystem. 

•  WinView/32  v2.4  is  used  for  camera  control  and  data  acquisition. 

•  IDL  5.5  and  ENVI  3.5  are  used  for  data  processing. 


Figure  5.4.  Pictures  of  LINUS.  On  the  left  side  is  the  electronics  section  with  the 
computer  the  amplifier  and  the  camera  pulser  and  cooler.  On  the  right  side  is  the  optics 
section,  which  also  includes  the  camera  body  and  the  scanner  mirror  assembly. 


3.  LINUS  Operation 

The  scanner  minor  is  set  to  a  desired  position  by  the  stepping  program.  Light 
coming  through  the  aperture  is  reflected  by  the  mirror  through  the  UV  filter  into  the 
optical  section  of  the  device.  The  filter  cuts  all  wavelengths  outside  its  response  region. 
The  radiation  is  analyzed  by  the  diffraction  grating  into  the  constituent  wavelengths  and 
is  recorded  by  the  camera.  The  camera  output  is  sent  to  the  host  computer  for  storage  and 
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subsequent  processing.  Once  image  acquisition  is  finished,  the  mirror  is  stepped  to  a  new 
position  and  the  process  is  repeated  until  the  whole  scene  is  imaged.  At  the  time  of  data 
acquisition  for  his  thesis,  mirror  stepping  and  image  acquisition  were  done  by  two 
separate  software  applications  (LabVIEW  and  WinView),  which  required  user  input. 
This,  in  addition  to  being  elaborate  and  time  consuming,  results  in  an  image  of  the  whole 
scene  that  is  extended  in  time  (as  much  as  20  min),  depending  on  the  camera  integration 
time  per  vertical  strip  of  the  scene.  Currently,  there  is  a  project  of  integrating  mirror 
scanning  and  data  acquisition  under  one  automated  process  using  LabVIEW. 

LINUS’s  product  is  a  dataset  consisting  of  three  components  for  each  pixel,  two 
spatial  coordinates  and  one  spectral  array.  The  vertical  spatial  coordinates  (y)  along  with 
the  spectral  components  for  each  horizontal  position  are  recorded  and  the  image  is 
scanned  horizontally  to  cover  the  scene  (see  figure  5.5).  The  data  can  be  thought  of  as  a 
hyperspectral  data  cube  shown  in  figure  5.6.  The  field  of  view  in  the  vertical  direction  (y 
axis)  is  0.5  degrees  and  is  recorded  in  512  pixels,  while  in  the  horizontal  direction  (x 
axis)  the  maximum  field  of  view  that  can  be  scanned  is  2  degrees,  or,  ±1  degree  from  the 
mirror  boresight  position.  The  horizontal  field  of  view  ( 0 ),  imaged  at  a  certain  mirror 
position,  depends  on  the  slit  width  setting  and  is  given  by  equation  5.1. 

0=2tan“1(^)  (5.1) 

where  0  is  in  radians,  w  is  the  slit  width  and  f=25cm  is  the  focal  length  of  the  primary 
objective  lens  of  the  optical  system  (see  figure  5.2  above).  The  number  of  wavelengths 
recorded  is  512  and  the  wavelength  resolution  depends  on  the  diffraction  grating  and  the 
slit  width  setting. 
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Figure  5.6.  Data  set  produced  by  LINUS  output. 
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VI.  EXPERIMENTAL  DATA  COLLECTION 


A.  GENERAL 

The  initial  schedule  to  field  test  LINUS  was  for  early  spring  of  2002.  Due  to 
component  manufacturing  and  integration  delays,  the  device  was  not  ready  until  late 
summer  of  the  same  year.  The  location  selected  for  the  first  data  acquisition  was  Lassen 
Volcanic  National  Park  with  the  main  site  being  the  Sulfur  Works  where  sulfuric 
emissions  occur.  The  Park  is  located  in  northern  California,  USA  as  shown  in  figure  6.1 
below.  The  Sulfur  Works  hydrothermal  area  (see  figure  6.2)  is  the  most  accessible  hot 
springs  area  in  Lassen  Park  and  it  is  thought  to  be  part  of  the  central  vent  system  of 
ancient  Mount  Tehama. 


LASSEN  VOLCANIC 
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Figure  6. 1 .  Map  of  Lassen  Volcanic  National  Park.  The  square  on  the  right  is 
magnification  of  the  area  included  in  the  left  box.  Locations  1,  2  and  3  are  the  spots 
where  the  measurements  took  place.  (After:  USGS  website,  2002). 
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Figure  6.2.  Some  of  the  Sulfur  Works  vents. 


B.  BRIEF  VOLCANIC  HISTORY  OF  LASSEN  PARK 

From  600,000  to  400,000  years  ago,  ancient  Mt.  Tehama  gradually  built  up 
through  countless  eruptions.  Later,  this  volcano  became  inactive  and  eroded  away  and  its 
main  vent  was  probably  what  is  now  the  park’s  Sulfur  Works.  Lassen  Peak  began  as  a 
volcanic  vent  on  Mt.  Tehama’s  northern  flank  and  was  formed  by  eruptions  about  27,000 
years  ago.  It  is  considered  the  world’s  largest  dome  volcano.  After  being  dormant  for 
27,000  years,  in  1914,  Lassen  Peak  started  a  seven-year  eruptive  activity.  The  most 
important  eruption  occurred  on  May  22,  1915  (see  figure  6.3)  (From:  USGS  website, 
2002). 
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Figure  6.3.  The  May  22,  1915  explosive  eruption  of  Lassen  Peak.  (From:  USGS 

website,  2002) 

C.  DATA  ACQUISITION 

1.  First  Trip 

The  first  trip  to  Lassen  took  place  between  21  and  23  August  2002.  Many  data 
sets  where  acquired  and  returned  back  to  NPS  for  evaluation  and  processing,  since  there 
was  no  capability  to  process  the  data  on  site  yet.  Evaluation  of  the  data  revealed  that  they 
were  essentially  useless,  and  further  investigation  showed  that  the  main  reason  was 
scanning -mirror  alignment  problems,  probably  caused  due  to  transportation  of  the  device. 
Nevertheless,  useful  experience  was  gained  regarding  the  transportation,  handling  and 
operation  of  LINUS  in  the  field. 

2.  Second  Trip 

The  second  trip  took  place  on  13  and  14  September  2002.  Evaluation  and 
processing  of  the  data  acquired  on  this  trip  forms  the  focus  of  this  thesis.  At  that  time,  the 
capability  of  scanning  across  the  scene,  although  available,  was  not  reliable.  Thus,  most 
of  the  images  taken  were  with  the  mirror  in  the  boresight  position. 

At  the  same  time  that  measurements  were  performed  with  LINUS,  other  data 
were  also  recorded  in  a  logbook  to  be  used  in  the  subsequent  processing  of  the 
measurements  and  modeling  with  MODTRAN.  Those  data  are: 

•  Date  and  time  of  day  (recorded  by  the  computer  when  saving  the  images). 
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•  Latitude,  longitude  and  elevation  given  by  a  hand-held  GPS  device. 

•  Inclination  with  respect  to  horizontal  and  heading  of  the  line-of-sight  of 
the  LINUS  optical  system. 

•  LINUS  optical  system  and  camera  settings:  slit  width,  camera  MCP 
voltage  and  integration  time. 

•  Turn  size  (in  degrees)  in  the  azimuthial  direction  when  the  controls  of  the 
tripod  that  supports  the  optical  system  were  used  to  perform  a  “mechanical 
scan”  and  step  size  for  the  few  cases  when  the  stepping  program  was  used 
to  perform  a  mirror- scan. 

•  Approximate  weather  conditions  (temperature,  wind,  clouds)  since  no 
instruments  were  available. 

•  Other  events  or  observations  considered  to  be  significant  for  the  follow  -on 
work  (e.g.  bad  images  for  specific  reasons,  etc  .). 

The  weather  conditions  were  ideal  for  taking  measurements  during  both  days, 
with  clear  skies,  mild  winds  and  an  average  temperature  of  about  22°C  (72°F).  The 
measurements  took  place  at  locations  1,  2  and  3  as  shown  on  the  map  in  figure  6.1.  Due 
to  heavy  road  construction  in  the  Park,  measurements  scheduled  to  take  place  at  two 
other  sites  were  not  possible.  The  geographical  coordinates  of  the  locations  are: 

•  Location  1:  40°26.92N,  121°32.14W,  elevation  6982  ft 

•  Location  2:  40°26.93N,  121°32.13W,  elevation  6993  ft 

•  Location  3:  40°26.88N,  121°31.98W,  elevation  8413  ft 

On  13  September,  LINUS  was  deployed  at  location  1  and  trial  images  were 
acquired  to  investigate  how  the  device  responded  to  combinations  of  various  inclinations 
and  azimuth  angles.  Slit  width  settings,  camera  MCP  voltage  settings  and  image 
integration  times  were  also  varied.  Most  of  the  measurements  were  taken  in  a  NE-E 
(northeast-east)  direction,  which  at  that  time,  was  about  90°  with  respect  to  the  sun’s 
position,  looking  over  the  area  where  SCL  was  emerging  (see  figure  6.2). 
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Figure  6.4.  Sulfur  Works  as  seen  from  location  1. 


On  14  September,  measurements  were  taken  at  locations  2  and  3.  At  location  2, 
LINUS  had  an  unobstructed  view  of  the  entire  Sulfur  Works  area  and  the  bulk  of  the  data 
were  acquired  there  (see  figure  6.3).  There  were  horizontal  scans  of  about  180°  in 
azimuth  with  the  Sulfur  Works  area  at  the  center  (±90°).  The  scans  were  performed 
manually  using  the  controls  of  the  tripod.  The  bulk  of  the  data  were  acquired  at  a  20° 
inclination  angle  with  respect  to  horizontal,  with  the  device  looking  at  the  sky  just  above 
the  tree-line. 
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Figure  6.5.  LINUS  setup  at  location  2. 


Location  3  was  at  the  parking  area  below  Lassen  Peak  away  from  Sulfur  Works. 
The  purpose  of  the  measurements  there  was  to  acquire  clear  sky  spectra.  These  will  be 
compared  to  those  from  locations  1  and  2  in  order  to  identify  expected  differences  in 
spectra  due  to  S02  presence  in  the  data  from  locations  1  and  2  with  no  or  much  less  S02 
at  location  3.  Images  with  the  input  aperture  covered  were  also  acquired  in  order  to 
record  the  dark  response  of  the  device. 

Difficulties  encountered  during  the  expeditions,  which  might  be  useful  to  mention 
for  the  next  expeditions,  are  as  follows: 

•  It  was  discovered  in  the  beginning  that  there  were  minor  light  leaks  into 
the  device  through  the  joints  between  the  housing  of  the  optics  section  and 
its  removable  lid,  which  provides  access  to  the  optics.  Although  minor, 
they  seemed  to  affect  the  data  under  certain  inclination  and  azimuthial 
angles  (position  of  the  optics  system  with  respect  to  the  sun).  The  light 
leaks  were  fixed  on-site  using  aluminum  tape  and  black  felt  where 
necessary. 

•  Under  certain  angles  of  the  optics  system  with  respect  to  the  sun,  light  was 
reflected  into  the  aperture  from  part  of  the  aluminum  frame  that  supports  it 
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on  the  tripod.  This  was  handled  by  covering  the  involved  portion  of  the 
frame  with  black  felt. 

The  device  is  quite  bulky  and  heavy  for  field  use  and  some  parts  of  it 
require  four  people  to  be  handled  properly  (loading-unloading  from  the 
vehicle  and  setup).  It  also  requires  a  paved  area,  so  that  the  transportation 
vehicle  can  access  it  and  the  device  can  be  deployed. 


D.  RAW  DATA 

During  the  two  days  of  LINUS  deployment,  625  data  files  were  acquired  taking 
up  303Mb  of  computer  memory  space.  Later,  these  data  were  evaluated  and  analyzed,  as 
described  in  chapter  VIII,  in  order  to  find  out  if  there  is  a  S02  signature  in  them  and 
whether  this  signature  can  provide  quantitative  estimates  for  the  gas  cone  entrations. 
Figure  6.6  shows  an  example  of  what  the  LINUS  operators  were  seeing  on  the  monitor  of 
the  device  at  the  time  of  the  data  files  acquisition.  Figure  6.7  is  the  same  as  figure  6.6 
with  the  difference  that  it  is  pseudocolored  in  order  to  get  a  feeling  for  the  radiation 
intensity  variations  along  the  spectrum.  Finally,  the  lower  part  of  figure  6.8  is  the  first 
step  in  the  data  file  processing,  which  gives  the  spectral  profile  of  the  recorded  image  by 
averaging  the  middle  256  lines  (from  128  to  384  strip  numbers).  The  units  are  Digital 
Numbers  in  all  figures. 
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Figure  6.6.  Screenshot  of  LINUS  monitor  at  the  time  of  data  file  acquisition.  This  is 
how  the  raw  data  looked  like  to  the  operators. 
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Figure  6.7.  The  previously  presented  data  file  in  pseudoc  olor,  which  gave  to  Linus 
operators  a  first  feeling  about  the  intensity  variations  along  the  spectrum  contained  in  the 

image. 
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VII.  MODELING  WITH  MODTRAN 


A.  GENERAL 

Apart  from  collecting  and  analyzing  experimental  data  using  LINUS,  it  was  also 
decided  to  simulate  the  data  with  a  radiative  transport  model  for  the  following  reasons: 

•  LINUS  had  never  been  used  in  the  field  before,  so  modeling  would  give 
an  indication  of  how  the  data  should  look  like  in  the  wavelength  -radiance 
space. 

•  Processing  of  the  modeled  data  would  provide  the  SO  2  presence  signature 
that  later  would  be  searched  for  in  the  experimental  data. 

•  Comparison  of  modeling  with  experimental  results  would  facilitate  S02 
detection  by  LINUS,  and  if  positive,  give  a  quantitative  estimate  of  the  gas 
concentration. 

A  sample  procedure,  which  can  be  applied  in  designing  the  modeling  of  electro  - 
optical  systems,  in  this  case  LINUS,  would  include  the  following  steps  (Anderson  et  al., 
1995): 

•  Selection  of  the  basic  problem  -  atmosphere  as  signature  or  contaminant. 
Atmosphere  is  modeled  as  a  signature  for  this  thesis. 

•  Selection  of  species,  viewing  platform,  necessary  geometric 

considerations  and  electro-optical  system  characteristics  to  be  employed. 
In  our  case,  the  species  of  interest  is  S02,  and  the  platform  is  LINUS 
operating  in  the  UV  with  a  specific  filter  response  (300  ±  10  nm)  and 
variable  viewing  geometries. 

•  Selection  of  the  modeling  algorithm  for  instrument  simulation  (line  by  line 
or  band  model),  based  on  ease  of  use  and  confidence  in  achieving  required 
accuracy.  The  algorithm  of  choice  for  this  case  was  the  band  model 
MODTRAN4. 

•  Selection  and/or  design  of  data  analysis  schemes  and  their 
implementation. 


B.  MODTRAN4 

In  the  early  1970’ s,  the  then  Air  Force  Cambridge  Research  Laboratory,  initiated 
a  program  to  develop  computer-based  atmospheric  radiative  transfer  algorithms.  The  first 
product  of  this  program  was  LOWTRAN  (LOW  resolution  TRANsmittance  code),  which 
was  released  in  1972.  This  limited  initial  effort  has  now  progressed  to  a  set  of  codes  and 
related  algorithms,  including  line-of-sight  spherical  geometry,  direct  and  scattered 
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radiance  and  irradiance,  non-local  thermodynamic  equilibrium  etc.,  that  contain 
thousands  of  coding  lines  in  hundreds  of  subroutines  with  improved  accuracy  and 
efficiency.  This  set  of  codes  includes,  among  others,  LOWTRAN,  FASCODE  (Fast 
Atmospheric  Signature  CODE)  and  MODTRAN  (MODerate  resolution  atmospheric 
TRANsmittance  and  radiance  code).  Each  code  relies  upon  a  catalogue  of  fundamental 
molecular  constants  as  contained  in  the  HITRAN  databases  (atlas  of  high -resolution 
molecular  spectroscopic  data)  (Anderson  et  al.,  1995). 

The  model  used  to  simulate  LINUS  experimental  data  is  MODTRAN4  (Berk  et 
al.,  1999).  MODTRAN4  is  currently  the  U.S.  Air  Force  standard  moderate  spectral 
resolution  model  for  atmospheric  transmission,  radiance  and  flux.  This  radiative  transport 
model  extends  from  the  thermal  infrared  through  the  visible  and  into  the  ultraviolet, 
covering  the  spectral  region  from  0.2pm  to  infinity  (50,000-Ocm1)  (Kneizys  et  al.,  1988). 
It  was  developed  jointly  by  Spectral  Sciences,  Inc.  and  the  Air  Force  Research 
Laboratory  /  Space  Vehicles  Directorate  (AFRL/VS).  MODTRAN  was  first  released  in 
1989  as  an  attempt  to  improve  the  spectral  resolution  of  LOWTRAN  7  (Kneizys  et  al., 
1988)  from  20cm'1  Full  Width  at  Half  Maximum  (FWHM)  with  5cm"1  calculation 
increments,  to  2cm  1  FWHM  with  1  cm  1  increments  (Berk  et  al.,  1989). 

MODTRAN  and  its  predecessors  have  been  used  extensively  over  the  last  quarter 
century  in  the  design  and  analysis  of  broadband,  multiband  and  short-wave  IR/Visible 
hyperspectral  imaging  sensors.  For  example,  remotely  sensed  spectral  imagery  of  the 
Earth’s  surface  can  be  used  to  its  fullest  advantage  only  when  the  influence  of  the 
atmosphere  has  been  removed  and  the  data  are  reduced  to  units  of  reflectance.  This  can 
be  achieved  with  MODTRAN  (Adle -Golden  et  al.,  1999). 

MODTRAN  uses  a  spherically  symmetric  atmosphere,  consisting  of 
homogeneous  layers  up  to  100km  altitude,  each  of  which  is  characterized  by  the  layer 
boundary  specification  of  temperature,  pressure  and  atmospheric  species  concentrations. 
It  uses  Snell’s  law  to  refract  a  line-of-sight  (LOS)  (Acharya  at  al.,  1999).  It  can  use  one  of 
the  six  available  geographical-seasonal  model  atmospheres  or  a  user  specified 
atmosphere.  The  atmospheric  species  whose  concentrations  can  be  user  -defined  are 
shown  in  the  table  7. 1  below. 
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No 

Species 

Other  name 

1 

water  vapor  (H20) 

2 

carbon  dioxide  (COO 

3 

ozone  (03) 

4 

nitrous  oxide  (N  ,0) 

5 

carbon  monoxide  (CO) 

6 

methane  (CH4) 

7 

oxygen  (07) 

8 

nitric  oxide  (NO) 

9 

sulfur  dioxide  (SO  2) 

10 

nitrogen  dioxide(NO  7) 

11 

ammonia  (NH3) 

12 

nitric  acid  (  HNO,) 

13 

CC13F 

CFC-11 

14 

CC19F9 

CFC-12 

15 

CC1F, 

CFC-13 

16 

cf4 

CFC-14 

17 

CHCIFt 

CFC-22 

18 

C,CTU 

CFC-113 

19 

C,C12F4 

CFC-114 

20 

c2cif5 

CFC-115 

21 

ciono2 

22 

hno4 

23 

CHChF 

24 

CC14 

25 

NA 

Table  7.1.  Molecular  species  for  user  defined  atmospheric  profiles.  SO2  is  number  9. 
(After  MODTRAN4  User’s  Manual,  1999) 


C.  SELECTION  OF  MODTRAN4  INPUT  CARDS  AND  VARIABLE  VALUES 

Although  MODTRAN  is  an  invaluable  scientific  tool,  for  the  novice  user,  it  is 
obscure  and  hard  to  use.  There  are  hundreds  of  variables  that  can  be  defined  by  the  user, 
which  are  only  very  briefly  described  in  the  user’s  manual  (Berk  et  al.,  1999). 
Furthermore,  the  user  interface  is  based  in  the  era  when  the  computers  were  initially 
developed  and  the  input  was  given  by  80  character  input  punch  cards.  Extensive  research 
in  the  existing  literature,  MODTRAN  documentation  and  published  applications,  did  not 
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reveal  many  details  about  most  of  the  variables  and  the  effects  they  have  on  the  output.  It 
seems  that  in  most  of  the  applications  slight  variations  of  the  standard  example  input  files 
(called  “cases”)  provided  with  the  MODTRAN  code  are  used.  For  the  purpose  of  this 
thesis,  none  of  the  sample  input  files  was  satisfactory  ev  en  with  extensive  modifications. 
Thus,  a  new  input  file  had  to  be  built.  In  order  to  understand  the  effects  of  certain 
variable  values  to  the  output  for  this  thesis,  extensive  trial  runs  had  to  be  performed. 
Some  of  the  results  are  presented  below,  with  the  purpose  of  documenting  what  was 
done,  and  helping  subsequent  users  avoid  repeating  the  same  procedure. 

1.  Selection  of  Input  Cards 

The  input  cards  used  for  this  modeling  case,  along  with  the  associated  variables 
and  their  values,  are  shown  in  appendix  A  Except  for  the  cards  required  by  the  model, 
there  are  a  number  of  optional  cards  that  are  used  to  specify  the  user  -defined  atmospheric 
profiles  and  geometry.  In  summary,  the  cards  used  are: 

•  Card  1 :  Main  radiation  transport  driver.  Required. 

•  Card  1A:  Radiative  transport  driver  (continued).  Required. 

•  Card  1 A2:  Spectral  data  and  sensor  response  function  file.  Optional. 

•  Card  2:  Main  aerosol  and  cloud  options.  Required. 

•  Card  2C:  User-defined  atmospheric  profiles.  Optional. 

•  Card  2C1:  User-defined  atmospheric  profiles.  Optional. 

•  Card  2C2:  User-defined  atmospheric  profiles.  Optional. 

•  Card  3STD:  Line-of-sight  geometry.  Required. 

•  Card  3A1:  Solar/Lunar  scattering  geometry.  Optional. 

•  Card  3A2:  Solar/Lunar  scattering  geometry.  Optional. 

•  Card  4:  Spectral  range  and  resolution.  Required. 

•  Card  5:  Repeat  run  options.  Required. 

2.  Selection  of  Variable  Values 

In  this  section,  the  values  selected  for  specific  variables  are  discussed.  Where  it 
was  clear  by  the  problem  what  the  values  should  be,  only  a  brief  justification  or 
description  is  given.  In  the  cases  where  it  was  deemed  that  further  investigation  was 
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needed  with  trial  runs,  all  the  other  variables  were  held  constant,  with  their  default  or  the 
previously  decided  values,  and  the  variable  in  question  was  changed.  For  some  of  these 
cases,  the  results  of  the  trial  runs  will  be  presented.  Examples  of  specific  input  files  used 
to  simulate  Lassen  measurements  are  presented  in  appendix  B. 


a.  Card  1 

MODEL=7:  It  specifies  that  user -defined  atmospheric  profiles  will  be 
provided  in  cards  2C,  2C1  and  2C2.  These  will  specify  the  way  the  atmosphere  will  be 
divided  in  layers  and  the  temperature  and  SCE  concentration  of  the  first  layer. 

IEMSCT=2:  The  code  has  to  execute  in  solar  radiance  mode  in  order  to 
model  the  radiance  reaching  the  sensor. 

IMULT=1:  The  spectral  region  of  interest  is  in  the  UV,  which  is 
dominated  by  Rayleigh  scattering.  Thus,  multiple  scattering  is  necessary.  Furthermore, 
the  solar  geometry  at  the  location  of  HI  (LINUS  latitude  and  longitude)  is  needed.  Trial 
runs  of  the  code  indicated  that  with  single  scattering,  radiances  are  systematically  smaller 
than  those  with  multiple  scattering,  with  the  maximum  difference  in  the  area  of  500nm 
where  the  radiance  curves  have  maxima.  The  effect  of  single  and  multiple  scattering  on 
the  radiance  for  an  extended  spectral  region  is  illustrated  in  figure  7.1  and  the  ratio  of 
single  scattering  divided  by  multiple  scattering  in  figure  7.2.  It  is  obvious  that  single 
scattering  gives  about  85%  less  radiance  than  multiple  scattering  in  the  spectral  region  of 
interest  (about  300nm). 
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Ratio  Radiance  (W/sr  err?  cm'1) 


Single  (IMULT=0)  and  Multiple  (IMULT=1)  Scattering  Plots 
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Figure  7.1.  MODTRAN  results  for  single  and  multiple  scattering. 


Ratio  of  Single  (IMULT=0t  to  Multiple  (IMULT=1)  Scattering 


Figure  7.2.  Ratio  of  single  divided  by  multiple  scattering  from  MODTRAN  results. 
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NOPRNT=0:  The  output  files  mainly  used  for  information  and  subs  equent 
processing  are  the  tape6  (TP6)  and  PLT.  If  NOPRNT=-l,  huge  TP8  files  are  created, 
which  generally  are  not  needed  but  require  a  lot  of  memory. 

b.  Card  1A 

DIS=T:  This  activates  the  DISORT  (Discrete  Ordinate  Radiative 
Transfer)  algorithm.  This  is  the  best  multiple  scattering  algorithm  implemented  in 
MODTRAN4,  which  together  with  the  following  DISAZM  and  NSTR  variables,  provide 
the  most  accurate  results  (Acharya  et  al.,  1999).  Figure  7.3  shows  that,  calculations 
without  using  the  DISORT  algorithm  deviate  from  the  better  calculations  involving 
DISORT,  particularly  for  wavelengths  smaller  than  292nm.  The  disadvantage  is  that  its 
use  increases  computational  time  by  a  factor  of  30  or  more. 


Ratio  of  Radiances  without  (DIS=R  and  with  (DIS=T)  DISORT  Algorithm 


Figure  7.3.  Effect  of  using  the  DISORT  algorithm  in  the  spectral  region  of  interest. 
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DISAZM=T:  Includes  azimuthial  dependence  in  the  line-of-sight  multiple 
scattering  of  solar  radiation.  It  is  needed  since  LINUS  can  scan  to  various  azimuthial 
directions  (headings). 

NSTR=8:  This  is  the  optimized  value  recommended  by  the  MODTRAN 
manual  for  most  applications.  There  are  no  specific  forward  scatterers  modeled  that 
would  require  the  value  of  16  except  from  the  default  ones. 

LSUN=F:  This  makes  use  of  the  default  solar  5cm  1  spectral  resolution 
irradiances.  Trial  runs  with  LSUN=T,  which  uses  the  1cm  1  binned  solar  irradiance  data 
and  ISUN=2,  which  assigns  FWHM=2cm  1  were  also  performed  to  investigate  how  this 
change  affects  the  output.  In  these  runs  it  was  set  LSUNFL=F  or  blank,  so  that  the 
Kurucz  solar  irradiance  database  was  used.  These  runs  showed  there  was  no  advantage  in 
using  LSUN=T  (see  figure  7.4),  especially  for  the  case  of  LINUS  which  has  a  spectral 
resolution  worse  than  9cm 1 . 


Ratio  of  LSUN=F  to  LSUN=T 


01280  0.285  0  290  0295  0.300  0.305  0  310  0315  0.320 
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Figure  7.4.  Solar  irradiance  comparison  for  5cm  1  (LSUN=F)  and  1cm 1  (LSUN=T) 
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spectral  resolution.  Results  for  Lassen  Park  with  MODTRAN  standard  atmosphere  2. 


C02MX=365.00000:  Use  the  recommended  CO2  mixing  ratio  (1999 
value)  instead  of  the  default  330  ppmv. 

LBMNAM=T:  In  order  to  read  the  newer  band  model  parameter  file 
BMP99_01.BIN  from  Card  1A2,  instead  of  the  default  BMP98_01.BIN. 

c.  Card  2 

IHAZE=  1 :  Sets  the  extinction  coefficient  to  rural  with  a  default  visibility 
of  23  km,  which  was  the  case  at  Lassen  Park  when  the  measurements  were  taken  . 

GNDALT=2. 13000:  This  is  the  altitude  in  km  at  Sulfur  Works  location  2 
in  Lassen  Park  (see  figure  6.1). 

d.  Card  2C 

ML=36:  Sets  the  number  of  layers  in  the  user  -defined  atmospheric  profile, 
which  in  this  case,  is  36. 

IRD1  =  1:  Read  data  for  user-defined  atmospheric  profiles  from  cards  2C1 

and  2C2. 


e.  Cards  2C1  and  2C2 

ZMDL  gives  the  altitude  of  layer  boundary  in  km.  After  several  trial  runs, 
it  was  found  that  MODTRAN,  whatever  the  initial  altitude  is  (say  x),  divides  the 
remaining  range  up  to  6km  in  6  layers  according  to  the  formula: 

GO-x 

6 

so  the  first  six  layers  (if  uniformly  thick)  are:  x,  x+y,  x+2y,  ...,  x+5y,  6.0.  After  6.0,  the 
layers  are  1km  thick  for  up  to  25km.  From  25  to  60km,  the  layers  are  5  km  thick,  from  60 
to  80km  they  are  10km  thick  and  from  80km  it  goes  directly  to  100km,  which  is  the 
boundary  of  the  last  layer.  For  the  Sulfur  Works  case,  the  first  layer  starts  at  2.13km.  By 
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the  geometry  of  the  Sulfur  Works  site  and  the  extent  of  the  plumes,  it  was  assumed  for 
modeling  purposes  that  the  plume  was  a  cylinder  with  a  diameter  of  approximately  250m 
(rising  vertically).  Furthermore,  by  assuming  that  LINUS  was  located  just  outside  the 
circumference  of  the  cylinder  and  looking  through  the  plume  at  an  average  inclination 
angle  with  respect  to  horizontal  of  20  degrees,  it  was  determined  that  the  first  layer, 
which  supposedly  contained  SO 2,  was  90m  thick.  Then,  the  boundaries  for  the  first  six 
layers  up  to  6km  were  calculated  to  be  2.130,  2.220,  2.976,  3.732,  4.488  and  5.244km. 
All  the  ZMDL  values  used  in  the  calculations  are  shown  in  sample  input  cards  presented 
in  appendix  B. 

T=2.200E+01:  This  is  the  temperature  in  Celsius  of  the  first  boundary.  For 
the  rest  it  was  set  to  T=0.000E+00  so  that  together,  with  subsequent  flags,  the  values  of 
the  MODTRAN  standard  model  atmosphere  2  (Mid -Latitude  Summer)  will  be  used. 

JCHAR(J)  (with  J=1  - 14)  =  2B22222222A222,  for  the  first  layer.  The  2 
means  use  the  default  values  for  the  standard  model  atmosphere  2.  JCHR(2)=B,  indicates 
ambient  temperature  in  degrees  Celsius.  JCHR(1 1)=A,  indicates  the  volume  mixing  ratio 
for  SO2  in  ppmv.  JCHAR(J)=22222222222222,  for  the  rest  of  the  layers,  which 
essentially  means,  that  except  for  the  first  layer  where  a  user  defined  temperature  and 
SO2  concentration  was  used,  for  the  res  t  of  the  layers,  default  values  from  the 
MODTRAN  standard  model  atmosphere  2  (Mid  -Latitude  Summer)  were  used. 

WMOL(J)=0.000E+00  for  J=  1  - 10  and  12.  Together  with  the  flags 
JCHAR(J)  makes  use  of  the  default  concentrations  from  the  standard  model  atmosphere  2 
of  the  gases  1- 10  and  12  shown  in  table  7.1  above. 

WMOL(l  1)  gives  the  concentration  of  the  SO2  in  ppmv.  Initial  runs  with 
the  regular  concentration  (WMOL(  1 1  )=  0.000E+00),  lOppmv  (WMOL(ll)=  1.000E+01), 
lOOppmv  and  l,000ppmv  were  performed  to  establish  a  SO 2  signature  in  the  spectrum. 
Subsequently,  runs  with  other  concentrations  were  performed  to  match  the  modeled  and 
experimental  spectral  profiles.  The  lOppmv  concentration  (0.721  atm -cm/km)  is  4.3xl04 
times  larger  that  that  of  the  MODTRAN  standard  atmosphere  2  under  the  same 
conditions.  The  SO2  concentration  as  a  function  of  altitude  is  given  in  figure  7.5. 
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SO 2  Concentration  Plot-  Results  from  MODTRAN  for  Standard  Atmosphere 2 


Altitude  (km) 

Figure  7.5.  SCT  concentration  as  a  function  of  altitude  for  the  MODTRAN  standard 

atmosphere  2. 


f.  Card  3 

HI  =2. 130:  This  is  the  initial  altitude  in  km  and  defines  the  position  of  the 

sensor. 

H2=blank  since  the  line-of-sight  is  a  slant  path  to  space. 

ANGLE=70:  This  is  90  degrees  minus  the  inclination  angle  of  the  line  -of- 
sight  of  the  instrument  with  respect  to  horizontal.  The  majority  of  the  measurements  were 
taken  at  an  inclination  of  20  degrees.  It  changes  depending  on  the  inclination  angle  where 
a  simulated  measurement  took  place. 

g.  Card  3A1 

IPARM=1:  Use  of  sensor  latitude  and  longitude  together  with  the  day  of 
the  year  and  time  of  the  day  to  specify  the  solar  geometry  with  respect  to  the  line-of-sight 
of  the  sensor. 
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IDAY=255:  This  is  the  Julian  day  that  the  measurements  took  place. 


h.  Card  3A2 

PARM  1=40.450:  This  is  the  latitude  of  location  2  at  Sulfur  Works.  It  is 
slightly  different  when  the  measurements  at  locations  1  and  3  are  modeled  (see  chapter 
VI). 

PARM2=121.530:  This  is  the  longitude  of  location  2  at  Sulfur  Works. 
Same  as  above  for  locations  1  and  3. 

TIME=  18.700:  This  is  the  Greenwich  time  of  day  when  a  modeled 
measurement  took  place.  It  changes  for  measurements  performed  at  different  times. 

PSIPO=290.000:  Shows  the  path  azimuth  or  heading  of  the  line-of-sight 
with  respect  to  the  North.  It  changes  depending  on  the  heading  where  a  simulated 
measurement  took  place. 

i.  Card  4 

VI  and  V2:  These  are  the  initial  and  final  frequencies  in  wavenumbers 
[cm1]  used  in  the  calculation  (spectral  range  of  the  simulation).  The  values  used  are  Vl  = 
20000  and  V2=50000  (0.2-0.5pm). 


D.  SIMULATION  NOTES 

1.  Variable  Values  Verification 

Due  to  the  structure  of  the  MODTRAN  input  file,  the  values  of  certain  variables 
earlier  in  the  file  restrict  or  define  the  use  or  not  of  subsequent  cards  and  the  values  of 
other  variables.  Thus,  caution  must  be  taken  so  that  the  values  used  are  compatible.  If 
not,  two  possibilities  normally  occur.  First,  the  code  does  not  run  and  gives  an  error 
message,  which  gives  the  user  an  indication  that  something  is  wrong.  Second,  the 
program  runs  without  giving  any  message  but  changes  the  incompatible  values  to 
compatible  ones  by  itself.  This  is  dangerous  and  can  produce  incorrect  results.  The  same 
effect  can  occur  due  to  a  misplaced  value  in  the  input  file,  which  is  possible  due  to  the 

strict  format  of  the  file.  For  example,  suppose  that  a  user  wants  to  run  MODTRAN  in  a 
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multiple  scattering  mode  using  the  DISORT  algorithm  and  azimuthial  dependence 
(DIS=T,  DISAZM=T,  NSTR=8).  If,  by  mistake,  IMULT=0  is  input  instead  of  the  correct 
1  in  the  previous  card,  the  program  will  assign  the  values  DIS=F,  DISAZM=F,  NSTR=0 
and  the  code  will  run  in  single  scattering  mode  without  giving  any  message.  A  good 
method  to  use  to  overcome  these  problems  is  to  check  the  variable  values  that  are  echoed 
in  the  tape6  output  file,  against  those  of  the  input  file  in  order  to  locate  any  possible 
changes. 

2.  Computational  Time 

The  use  of  the  DISORT  multiscatering  algorithm  with  azimuthial  dependence,  in 
order  to  produce  more  accurate  results,  increases  the  computational  time  greatly.  Runs 
without  using  DISORT  require  2-4  minutes  on  a  personal  computer  with  an  AMD 
850MHz  CPU,  depending  on  the  values  of  some  variables.  The  same  runs  with  the 
DISORT  algorithm  require  about  1.5-2  hours,  which  is  about  30  times  longer  run  times 
on  average.  This  can  be  a  problem  if  a  large  number  of  simulations  is  required  in  a  short 
period  of  time. 


3.  Spectral  Resolution 

In  the  initial  trial  runs  where  SO2  absorption  features  were  investigated, 
MODTRAN  was  run  at  the  finest  possible  resolution  of  1cm  1  bins  with  2cm"1  FWHM. 
During  those  runs  it  was  found  that  if  the  spectral  lower  limit  VI  was  greater  than 
22,681cm  1  (lower  that  0.441pm),  the  program  did  not  run  with  spectral  bins  of  lcm"1. 
Instead,  it  required  resolution  equal  or  greater  than  5cm  1  and  ran  in  LOWTRAN  mode. 
For  this  reason,  the  spectral  region  used  in  the  initial  simulations  was  50,000 -20,000cm  1 
(0.2-0. 5pm).  The  output  files  were  later  reduced  to  the  spectral  region  of  interest  (280- 
320nm).  Due  to  the  inverse  relation  between  spatial  frequency  (cm1)  and  wavelength,  the 
spectral  resolution  over  extended  spectral  regions  is  not  constant  but  gets  finer  as  the 
wavelength  decreases.  For  the  300nm  region,  the  resolution  is  0.009nm  if  the  spectral 
bins  are  lcm1. 
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The  actual  resolution  of  LINUS  depends  on  the  slit  width  and  the  exact  values 
were  not  known  at  the  time  this  thesis  was  written.  Estimates  indicate  that  it  should  be 
from  0.4  to  lnm.  When  LINUS  measurements  were  simulated,  runs  at  various  resolutions 
were  performed  in  order  to  reproduce  profiles  with  similar  spectral  details  as  the 
measured  ones.  The  finest  resolution  used  was  9cm  1  (0.08nm  at  300nm)  and  was  based 
on  the  assumption  that  the  512  pixels  of  the  camera  focal  plane  sample  one  wavelength 
per  pixel  in  the  region  278-3 18nm  where  LINUS  operates  (see  wavelength  calibration  in 
chapter  VIII).  The  results  were  much  more  detailed  than  the  measured  spectral  profiles. 
After  several  trials  it  was  found  that  simulations  with  55cm  1  bins  (0.5nm  at  300nm) 
reproduced  results  similar  in  detail  with  the  LINUS  measurements. 
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VIII.  DATA  ANALYSIS 


A.  SIMULATED  DATA 

1.  LINUS  UV  Filter 

As  mentioned  in  chapter  V,  after  the  radiation  is  reflected  on  the  scanning  mirror, 
it  enters  the  optics  section  of  the  device  through  a  filter  with  the  response  given  by  figure 
5.3.  The  vendor  of  the  filter,  Omega  Optical,  Inc.,  was  contacted  in  order  to  provide  the 
filter  data  in  electronic  form  so  they  could  be  used  in  the  simulation  but  they  did  not 
respond  to  the  request.  The  filter  response  curve  was  then  digitized  and  fitted  by  the 
Gaussian  given  by  equation  8.1,  in  units  of  nm,  and  it  is  plotted  in  figure  8.1.  This  fit  was 
then  used  to  simulate  the  effect  of  the  filter  on  the  MODTRAN  output  data. 

,,r-298.43,2/o 

f(x)  =  15.11  e  448  (8.1) 


280  290  300  310 

Wavelength  (nm) 

Figure  8.1.  The  continuous  line  shows  the  Gaussian  fit  of  the  LINUS  filter  response. 
Diamonds  show  data  points  taken  from  the  real  filter  response  curve  presented  in  figure 

5.3. 
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2.  S02  Presence  Effects  and  Spectral  Signature 

MODTRAN  was  run  in  order  to  find  out  how  the  presence  of  S02  in  the  plume 
affects  the  spectral  profile  and  to  establish  a  signature  that  would  be  searched  for  in  the 
experimental  data.  Runs  for  normal  atmosphere,  and  atmospheres  with  lOppmv, 
lOOppmv  and  lOOOppmv  of  S02  were  performed.  Since  at  the  time  this  thesis  was 
written,  the  spectral  resolution  of  LINUS  was  not  known  precis  ely,  the  simulations  were 
done  in  two  spectral  resolutions.  The  first  is  0.08nm  (9cm  1  at  300nm)  and  represents  the 
best  achievable  resolution  assuming  that  each  of  the  512  camera  pixels  records  a  different 
spectral  channel  (no  over-  or  under-sampling)  in  the  range  278.15-319.20nm,  where 
LINUS  operates  as  measured  during  calibration  with  platinum  spectra.  The  second  is 
0.5nm  (55cm  1  at  300nm),  which  is  the  estimated  device  resolution  taking  into  account 
the  device  slit  width  and  optics  limitations.  The  results  presented  in  figures  8.2  through 
8.8  are  for  the  Sulfur  Works  location  2,  with  a  heading  of  290  degrees,  an  inclination  to 
horizontal  of  20  degrees  at  10:45  a.m.  The  lines  in  the  plots  are  fitted  curves  to  the  data 
points,  which  are  plotted  with  dots. 


a.  Figures  8.2  through  8.4 

Figures  8.2  through  8.4  show  the  simulated  spectral  profiles  for  the 
various  S02  concentrations  that  reach  the  aperture  of  LINUS  before  they  pass  through  the 
filter.  Figures  8.2  and  8.3  are  for  the  same  resolution  of  0.08nm,  but  the  second  picture 
has  a  smaller  scale  in  the  radiance  axis  in  order  to  show  the  details  at  lower  wavelengths. 
Figure  8.4  shows  the  profiles  for  0.5nm  resolution.  The  observations  deduced  are  as 
follows: 

•  The  radiation  is  rapidly  reduced  as  the  wavelength  is  reduced  and  there  is 
practically  no  signal  below  295nm  even  for  the  regular  atmosphere. 

•  The  profile  of  the  low  S02  concentration  of  lOppmv  is  similar  to  that  of 
the  regular  atmosphere  with  minor  differences.  The  profile  of  the 
lOOppmv  concentration  is  significantly  different  than  that  of  the  regular 
atmosphere  and  strong  S02  absorption  features  can  be  identified.  The 
lOOOppmv  concentration  diminishes  the  signal  so  much  that  it  makes  it 
practically  zero,  so  that  presence  of  concentration  of  this  order  in  the  line- 
of-sight  of  the  device  would  be  immediately  evident. 
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Radiance  (W/cnr  sr  cm 


•  There  is  a  significant  difference  in  the  structure  of  the  profiles  with 
resolution.  The  0.08nm  profiles  have  many  more  details  than  that  of 
0.5nm,  facilitating  the  identification  of  absorption  characteristics.  In  the 
0.5nm  resolution,  the  details  are  smoothed  out  and  in  order  to  identify 
differences  in  the  profiles,  the  difference  in  concentrations  has  to  increase 
significantly. 


VARIOUS  S02  CONCENTRATIONS  WITHOUT  FILTER  EFFECT  -  RESOL.  0.08nm 
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Figure  8.2.  Simulated  sky  radiance  spectral  profiles  for  various  SCL  concentrations  at 
Sulfur  Works  without  LINUS  filter  effect.  Resolution  0.08nm. 
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VARIOUS  SQ2  CONCENTRATIONS  WITHOUT  FILTER  EFFECT 

RESOL.  0.08nm  -  SMALL  SCALE 


290  295  300  305  310  315  320 

Wavelength  (nm) 


Figure  8.3.  Simulated  sky  radiance  spectral  profiles  for  various  SO  2  concentrations  at 
Sulfur  Works,  without  LINUS  filter  effect,  at  smaller  scale  than  figure  8.2.  Resolution 

0.08nm. 


VARIOUS  SQ2  CONCENTRATIONS  WITHOUT  FILTER  EFFECT  -  RESOL.  0.5nm 


Figure  8.4.  Simulated  sky  radiance  spectral  profiles  for  various  SO  2  concentrations  at 
Sulfur  Works  without  LINUS  filter  effect.  Resolution  0.5nm. 
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b.  Figures  8.5  and  8.6 

Figures  8.5  and  8.6  show  the  simulated  spectral  profiles  for  the  var  ious 
SO2  concentrations  after  they  have  past  through  the  LINUS  filter.  Essentially,  these 
figures  are  produced  from  the  data  used  to  plot  figures  8.2  and  8.4  after  applying  the 
Gaussian  effect  shown  in  figure  8.1.  The  observations  deduced  from  these  figures  are: 

•  It  is  obvious  that  the  filter  suppresses  the  spectral  profile  above  312nm 
where  significant  SO2  absorption  effects  are  observed  in  the  unfiltered 
data.  There  is  practically  no  effect  below  295nm  from  the  filter,  since  the 
available  radiation  is  already  diminished  by  itself. 

•  The  signal  passing  through  the  lOOOppmv  concentration,  already  small 
without  the  filter  effect,  after  the  filter  is  so  small  that  it  does  not  show  up 
with  the  existing  scale  of  the  radiance  axis.  For  this  reason,  it  was  not 
plotted  in  this  set  of  figures. 

•  The  filter  modulates  the  profiles  and  gives  them  a  Gaussian -like  shape. 
Although  the  filter  has  maximum  transmittance  at  298.4nm,  the  maximum 
of  the  spectral  profiles  is  at  about  305nm  due  to  the  rapid  increase  of  the 
signal  as  the  wavelength  increases. 

•  As  with  the  previous  figures,  the  0.08nm  resolution  profiles  have 
significantly  more  detailed  structure  than  those  of  the  0.5nm.  With  the 
latter  resolution,  much  of  the  detail  is  lost  due  to  under -sampling  of  the 
spectral  region. 

•  There  is  a  number  of  well-defined  S02  absorption  characteristics  in  the 
spectral  profiles  of  the  S02  concentrations,  with  much  more  profound 
those  of  the  lOOppmv.  These  characteristics  are  greater  in  number  and 
magnitude  in  the  0.08nm  resolution  plot  and  are  diminished  in  the  0.5nm 
resolution  plot.  In  the  first  plot  (figure  8.5)  they  are  observed  at  297.9, 
299.9,  300.6,  302.4,  302.9,  304.0,  304.6,  306.3,  308.6  and  310.7nm.  In  the 
second  plot  (figure  8.6)  they  appear  at  298.0,  299.9,  308.6  and  310.7nm. 
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Figure  8.5.  Simulated  spectral  profiles  for  various  SO2  concentrations  at  Sulfur  Works 
with  LINUS  filter  effect.  Resolution  0.08nm 
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Figure  8.6.  Simulated  spectral  profiles  for  various  S02  concentrations  at  Sulfur  Works 
with  LINUS  filter  effect.  Resolution  0.5nm 
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c.  Figures  8.7  and  8.8 

Figures  8.7  and  8.8  show  the  ratios  of  the  SO2  containing  spectral  profiles 
divided  by  the  profile  of  the  regular  atmosphere  in  order  to  find  out  the  SO  2  presence 
signature.  The  observations  from  these  figures  are: 

•  Both  0.08  and  0.5nm  resolutions  give  similar  results,  but  the  absorption 
characteristics  in  the  former  are  well  defined  by  many  data  points,  while  in 
the  latter,  they  are  barely  defined  by  3-4  data  points.  This  indicates  that  the 
resolution  of  0.5nm  is  marginal  and  lo  wer  resolutions  (larger  than  0.5nm) 
are  not  able  to  identify  SO2  absorption  characteristics. 

•  The  troughs  are  due  to  the  SO2  absorption  bands  with  the  stronger  effect 
appearing  at  299. 9nm  for  the  lOppmv  concentration  and  299.9,  304.4  and 
306. 5nm  for  the  lOOppmv.  Furthermore,  the  total  signal  is  diminished 
especially  at  smaller  wavelengths.  The  greater  the  concentration,  the 
smaller  the  overall  signal.  At  lOOOppmv,  there  is  practically  no  signal  as 
was  also  observed  in  the  previous  plots  of  the  radiance  space. 

•  The  absorption  characteristics  are  in  agreement  with  those  observed  in  the 
previous  figures  of  the  wavelength-radiance  plots. 

•  The  absorption  characteristics  revealed  by  the  ratios  are  in  a  very  good 
agreement  with  the  SO2  absorption  cross-section  measurements  presented 
in  chapter  III,  figures  3.2  and  3.7. 

•  The  ratio  is  very  sensitive  to  S02- induced  effects  in  the  spectral  profiles. 
Even  for  the  low  concentration  of  lOppmv,  which  does  not  induce 
significant  changes  with  respect  to  the  profile  of  the  regula  r  atmosphere  in 
the  radiance  space,  the  ratio  clearly  reveals  its  presence.  This  indicates  that 
the  ratio  with  clear  sky  spectra  is  a  valuable  tool  giving  a  definite  S02 
signature  even  for  small  concentrations. 
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Figure  8.7.  Ratios  of  SO2  concentrations  to  regular  atmosphere  at  0.08nm  resolution. 
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Figure  8.8.  Ratios  of  SO?  concentrations  to  regular  atmosphere  at  0.5nm  resolution. 
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3.  Spectral  Profile  Structure  Investigation 

In  an  effort  to  identify  the  origins  of  the  observed  spectral  profile  structur  e  in  the 
region  of  interest,  additional  MODTRAN  simulations  were  performed.  In  these 
simulations,  the  effects  of  the  oxygen  (CF)  and  ozone  (O3)  molecules  were  investigated 
since  these  are  considered  to  be  the  most  significant  molecules  in  the  UV  region. 
Comparisons  of  the  results  of  the  MODTRAN  standard  atmosphere  2  (Mid  -Latitude 
Summer)  to  those  of  the  same  atmosphere  that  contained  no  O2,  no  O3  and  no  O2  and  O3 
at  the  same  time  were  performed.  Solar  irradiance  measurements  at  the  top  of  the 
atmosphere  (TOA)  were  also  examined.  The  following  conclusions  were  reached  from 
those  simulations: 

•  Oxygen  has  no  effect  in  the  structure  of  the  spectral  profile  in  the  region  of 
interest. 

•  Ozone  greatly  diminishes  the  radiance  by  a  varying  way  as  illustrated  in 
figures  8.9  through  8.11  below.  In  figure  8.9,  the  radiances  of  regular 
atmosphere  and  atmosphere  without  O3  are  plotted.  For  illustration 
purposes,  a  moving  average  is  also  plotted.  In  figure  8.10,  the  ratio  of  the 
regular  atmosphere  to  atmosphere  with  no  Q  in  bgarithmic  scale  is 
presented  in  order  to  illustrate  the  varying  effect  of  O3  with  wavelength.  In 
figure  8.11,  the  same  ratio  in  linear  scale  is  plotted  in  order  to  highlight  the 
region  above  300nm.  It  is  evident  that  at  smaller  wavelengths,  O3  does  not 
induce  a  structure  to  the  spectral  profile  up  to  313nm  except  that  it 
significantly  and  monotonically  decreases  the  signal. 


SECTRUM  OF  REGULAR  ATMOSPHERE  AND  ATMOSPHERE  WITH  NO  0, 


Spectrum  of  regular  atmosphere  and  atmosphere  with  no  O3. 
11 


Figure  8.9. 
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RATIO  OF  REGULAR  AT  MQSPHER  E  T  0  ATMOSPHERE  WfTH  NQQ, 
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Figure  8.10.  Ratio  of  radiances  in  log -scale  of  regular  atmosphere  to  atmosphere  with 

no  03. 
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Figure  8.1 1.  Linear  ratio  of  radiances  of  regular  atmosphere  to  atmosphere  with  no  O3. 
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Radiance  of  Reg.  Atm.  (W/cm2  sr  cm-1) 


In  figure  8.12,  the  solar  irradiance  out  of  the  atmosphere,  as  given  in  the 
corrected  Kurucz  database  (Berk  et  al„  1999),  together  with  the  radiance 
of  regular  atmosphere  at  2.13km  altitude,  are  plotted.  Due  to  the 
significant  scale  change  that  the  atmospheric  radiance  experiences  in  this 
region,  only  a  small  spectral  interval  is  plotted.  For  illustration  purposes,  a 
moving  average  w  as  used.  There  is  good  agreement  in  the  overall  structure 
of  the  two  curves.  Thus,  the  main  cause  for  the  observed  spectral  structure 
in  the  spectral  region  of  interest  is  the  solar  irradiance  before  it  enters  the 
atmosphere. 


REGULAR  ATMOSPHERE  RADIANCE  AND  SOLAR  IRRADIANCE  COMPARISON 
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Figure  8.12.  Spectrums  of  solar  irradiance  out  of  the  atmosphere  and  regular 
atmosphere  radiance  at  2.13km  altitude. 


B.  EXPERIMENTAL  DATA 

1.  Wavelength  and  S02  Calibration  of  LINUS 

Soon  after  the  acquisition  of  the  experimental  data,  a  wavelength  calibration  was 
performed  in  the  lab  in  order  to  assign  a  wavelength  scale  to  the  raw  data.  The  calibration 
was  performed  using  a  platinum  lamp  and  the  results  are  shown  in  figure  8.13. 
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WAVELENGTH  CALIBRATION  OF  LINUS 


Figure  8.13.  Results  of  wavelength  calibration  of  LINUS  using  a  platinum  lamp. 


An  SO2  calibration  was  also  performed  using  a  photocell  with  various  pressures 
of  SO2  and  a  deuterium  lamp  as  the  UV  source.  The  results  with  a  slit  width  of  0.1 1mm 
are  presented  in  figure  8.14.  The  S  (T  absorption  characteristics,  the  reduction  in 
transmission  and  the  fine  structure  of  the  curves,  are  evident.  It  should  be  noted  that  the 
fine  structure  did  not  appear  in  the  measurements  performed  with  a  0.55mm  slit  width, 
due  to  the  reduced  resolution  of  the  device  at  that  slit  width.  Each  plotted  curve  is  the 
ratio  of  the  radiance  recorded  by  the  device  at  the  specific  SO2  pressure,  divided  by  a 
subsequent  measurement  with  the  photocell  completely  evacuated.  For  more  details  on 
the  two  calibrations,  see  Gray,  2002. 
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2.  Data  Analysis 

For  each  collected  image  a  wavelength  assignment  was  performed  according  to 
the  results  of  the  calibration  mentioned  earlier.  Then,  a  spectral  profile  was  extracted  by 
averaging  the  radiance  values  of  the  middle  256  pixels  along  the  vertical  axis  of  the 
image  (spatial  y  axis)  as  explained  and  illustrated  in  chapter  VI,  figure  6.8.  Subsequently, 
there  were  two  approaches  in  the  data  analysis:  an  individual  image  analysis  and  an 
ensemble  analysis. 


a.  Individual  Image  Analysis 

Initially,  the  selection  of  specific  images  was  done  according  to  several 
criteria,  such  as  the  quality  of  the  image,  no  over  or  under  exposure,  no  objects  in  the 
image,  strong  spectral  profile  etc.  From  those  images,  those  taken  with  the  line-of-sight 
over  the  gas  vents  and  those  of  the  clear  sky  taken  at  Lassen  Peak  were  finally  selected. 
Images  taken  with  the  aperture  covered  were  also  selected  in  order  to  evaluate  the 
response  of  the  device  with  no  incoming  signal. 

The  plots  of  the  spectral  profiles  of  the  selected  images  are  presented  in 
appendix  C.  Figure  C.l  shows  the  spectral  profiles  of  all  selected  images.  Figures  C.2  and 
C.3  show  spectral  profiles  of  images  taken  over  the  vents  at  Sulfur  Works  from  location  2 
and  figure  C.4  shows  those  taken  at  Lassen  Peak  (clear  sky  spectra).  In  figure  C.5,  the 
data  were  acquired  with  the  aperture  of  LINUS  covered.  In  figures  C.2  through  C.5,  the 
red  curve  represents  the  average  of  all  the  plotted  curves  in  each  figure. 

A  plot  of  all  the  averages  is  shown  in  figure  8.15  below.  The  differ  ences  in 
radiance  are  due  to  different  slit  widths,  acquisition  times,  inclinations,  headings  and  time 
of  the  day  the  images  were  acquired.  The  sw_012  through  sw_020  files  are  the  ones  with 
the  aperture  covered  (see  also  figure  C.5).  It  is  obvious  that  their  spectral  profile  is 
practically  flat  and  does  not  induce  any  structure  to  the  signal.  For  that  reason  and  since 
ratios  were  mainly  used  in  the  subsequent  analysis,  no  action  to  eliminate  the  device 
“dark”  response  was  taken. 
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AVERAGES  OF  THE  SELECTED  DATA  GROUPS 


Figure  8.15.  Plot  of  the  averages  from  figures  C.2  through  C.5. 


Figure  8.16  shows  a  comparison  in  radiance  space  of  selected  spectral 
profiles  with  strong  structure  taken  over  the  Sulfur  Works  vents,  where  the  presence  of 
SO2  was  investigated,  divided  by  a  Lassen  Peak  profile,  which  is  clear  atmosphere  data. 
The  most  notable  differences  in  structure  appear  at  297.7  and  299.7nm,  where  two  peaks 
of  the  Lassen  Peak  profile  do  not  appear  in  the  Sulfur  Works  profile.  These  two 
characteristics  are  in  agreement  within  0.3nm  with  the  SO2  absorption  cross-section 
characteristics  measured  at  298.0  and  300.0nm  (chapter  III,  figure  3.3)  and  within  0.2nm 
with  the  simulation  results  at  297.9  and  299. 9nm  (figure  8.5  above). 
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EXPERIMENTAL  DATA  AT  SULFUR  WORKS  SW  OOP.  D80.2D6 

AND  LASSEN  PEAK  SW  503 


Figure  8.16.  Comparison  of  spectral  profiles  with  strong  structure  over  the  Sulfur 
Works  vents  (possible  S02  presence)  with  Lassen  Peak  (clear  atmosphere). 


Furthermore,  the  profiles  with  the  most  significant  spectral  structure  from 
Sulfur  Works  were  divided  by  one  from  Lassen  Peak,  in  an  effort  to  find  out  whether  a 
signature  of  SCb,  as  specified  by  the  simulations,  was  present.  The  results  are  presented 
in  figure  8.17  below.  The  observations  from  this  figure  are  as  follows: 

•  The  overall  shape  of  the  curves  has  similarities  to  those  of  the  SO  > 
calibration  (figure  8.14),  the  difference  being  that  the  calibration  curves 
are  shifted  about  13nm  to  the  left.  In  the  calibration  case,  the  radiation 
source  was  a  deuterium  lamp  while  at  Lassen  Park,  it  was  the  sun. 

•  There  is  a  structure  to  the  curves  with  troughs  at  297.8,  299.8,  300.6, 
301.6,  303.7,  305.8,  308.4  and  311.1nm.  From  those,  the  troughs  at  301.6, 
303.7  and  305. 8nm  are  very  profound  while  the  rest  are  weak.  The 
troughs,  which  agree  with  the  S02  absorption  characteristics  specified  by 
the  MODTRAN  simulations,  are  at  297.8,  299.8,  300.6  and  308.4nm.  The 
rest,  including  the  most  profound,  do  not  agree  with  those  that  were 
simulated.  The  differences  are  due  to  a  right  shift  of  about  0.5-0.6nm  of 
the  simulated  data. 
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RATIOS  OF  SELECTED  SULFUR  WORKS  TO  LASSEN  PEAK  DATA 


280  285  290  295  300  305  310  315  320 

Wavelength  (nm) 


Figure  8.17.  Ratios  of  strong  Sulfur  Works  spectra  (possible  SOi  existence)  with  strong 
Lassen  Peak  spectra  (clear  atmosphere). 


The  ratios  of  the  averaged  spectral  profiles  are  plotted  in  figure  8.18 
below.  The  detail  of  the  spectral  structure  of  the  previous  figure  is  diminished  but  not 
lost.  Ratios  of  profiles  with  a  less  strong  structure  are  also  plotted  in  figure  8.19  for 
comparison.  These  ratios  have  no  detail  in  their  structure  compared  to  figures  8.17  and 
8.18. 
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RATIOS  OF  AVERAGED  DATA  FROM  SULFUR  WORKS  LOC.  2  TO  LASSEN  PEAK 


Figure  8.18.  Ratios  of  the  Sulfur  Works  averages  (blue,  green)  to  the  Lassen  Peak 
average  (red)  show  n  in  figure  8.15  above. 


RATIOS  OF  SELECTED  SULFUR  WORKS  TO  LASSEN  PEAK  DATA 


Figure  8.19.  Ratios  of  Sulfur  Works  (possible  SO2  existence)  with  Lassen  Peak  spectra 
(clear  atmosphere).  Spectral  profiles  that  do  not  have  strong  structure  were  used. 
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b.  Ensemble  Analysis 

In  this  approach,  the  data  were  treated  as  a  whole  in  an  effort  to  find  bold 
characteristics  of  SO2  absorption.  One  of  the  products  of  this  effort  is  illustrated  in  figure 
8.20.  The  vertical  axis  shows  all  the  images  collected  at  Lassen  Park  represented  as  lines. 
In  order  to  produce  a  line,  the  middle  256  spatial  lines  of  each  image  were  averaged  as 
described  earlier.  Horizontal  black  lines  represent  images  that  were  rejected  as  “bad”,  in 
order  not  to  induce  undesirable  effects  to  the  statistics  of  the  whole  image.  The  average 
spectrum  of  all  images,  along  with  a  simulated  spectrum  of  clear  atmosphere,  are  also 
plotted.  Finally,  the  figure  includes  one  of  the  results  of  the  SO 2  calibration.  Observations 
similar  to  those  deduced  from  the  individual  image  analysis  described  earlier  can  be 
identified.  The  spectral  character  of  all  images  is  evident.  No  bold  SO  2  absorption 
characteristics  could  be  identified  by  this  approach. 
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Figure  8.20.  Summary  of  the  experimental  measurements  with  a  simulated  clear 
atmosphere  and  SO 2  calibration  spectral  profiles . 
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C.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  SIMULATED  DATA 

In  an  effort  to  evaluate  how  well  simulated  and  experimental  data  agree  with  each 
other,  and  if  SO2  presence  characteristics  deduced  by  the  simulations  exist  in  the 
experimental  data,  many  comparisons  were  performed.  The  comparisons,  as  above, 
included  selected  images  against  selected  simulated  data  and  comparisons  of  all  the  data 
as  a  whole  as  described  in  the  ensemble  approach  above.  Some  of  these  comparisons  in 
the  radiance  space  are  presented  in  figures  8.21  through  8.26. 

The  general  observations  from  these  graphs  are: 

•  There  is  a  good  agreement  between  simulated  and  experimental  data  from 
about  290  to  306nm.  From  306nm  and  above,  although  the  profile 
structures  are  similar,  the  experimental  profile  becomes  wider  and  the 
peaks  and  troughs  are  in  different  wavelength  positions.  This  is  attributed 
to  effects  of  the  filter  of  the  device  to  the  recorded  radiation,  which  could 
not  be  reproduced  by  the  filter  simulation.  As  illustrated  in  figure  5.3,  the 
filter  response  curve  is  slightly  tilted  towards  the  smaller  wavelengths  and 
the  Gaussian  fit  did  not  reproduce  this  effect  (see  figure  8.1). 

•  The  shape  of  the  experimental  curves  along  the  radiance  axis  depends 
upon  the  acquisition  time  and  slit  width,  while  the  simulated  curves  do 
not. 

•  The  simulated  curves  have  zero  radiance  outside  the  response  region  of 
the  simulated  filter,  while  the  experimental  curves  have  non -zero  values  in 
this  region  due  to  the  dark  response  of  the  device  and  possibly,  stray 
radiation  coming  into  the  device  from  paths  other  than  its  line-of-sight. 

Figures  8.21  and  8.22  show  the  simulated  and  measured  spectral  profiles  of  clear 
atmosphere  (without  SCF).  The  simulated  data  are  in  the  two  resolutions  used  in  the 
simulations.  It  is  evident  that  the  details  of  the  experimental  spectral  profile  match  those 
of  the  0.5nm  resolution  much  better  that  those  of  the  finer  0.08nm  resolution.  The 
conclusion  from  these  graphs  is  that  LINUS  with  a  0.11mm  slit  width  operates  at  about 
0.5nm  resolution. 
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EXPERIMENTAL  VS  SIMULATED  (AT0.08nm  RES.)  CLEAR  ATMOSPHERE  SPECTRA 


Figure  8.21.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere, 
simulated  data  are  in  0.08nm  resolution. 


EXPERIMENTAL  VS  SIMULATED  (AT0.5nm  RES.I  CLEAR  ATMOSPHERE  SPECTRA 


Figure  8.22.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere.  The 
simulated  data  are  in  0.5nm  resolution. 


89 


Radiance  (DN)  SW_503  Radiance  (DN)  -  SW  503 


Figure  8.23  shows  two  spectral  profiles  looking  over  the  vents  at  Sulfur  Works, 
compared  to  regular  atmosphere  and  lOppmv  S02  simulated  spectral  profiles.  The  most 
noticeable  difference  in  the  simulated  data  at  300nm  agrees  with  the  behavior  of  the 
experimental  curves  in  that  area. 


SULFUR  WORKS  VS  SIMULATED  (AT0.5nm  RES.)  SPECTRA 


2800 

2600 

LO 

s, 

2400  § 

O 
CO 

2200  gl 

j/j 

2000  z 

Q 

0) 

o 

1800  S 

B 
cc 

1600 


Figure  8.23.  Measured  and  simulated  spectral  profiles  of  clear  atmosphere  and 
simulated  with  lOppmv  SCF  atmosphere.  The  simulated  data  are  in  0.5nm  resolution. 


Furthermore,  ratios  between  experimental  and  simulated  data  were  investigated. 
The  simulated  data  were  remapped  to  the  experimental  data  wavelength  coordinates 
using  interpolation,  so  that  the  ratios  could  be  calculated.  Some  of  the  results  are 
presented  in  figures  8.24  through  8.26. 

Figure  8.24  shows  the  ratios  of  two  images  taken  at  Sulfur  Works  (SW_080, 
SW_206)  and  one  at  Lassen  Peak  (SW_503)  divided  by  simulated  clear  atmosphere 
spectra.  The  first  two  images  are  taken  over  the  vents  while  the  third  is  a  clear 
atmosphere  image.  The  structure  of  all  the  curves  looks  nearly  id  entical  with  the  most 
noticeable  difference  at  300nm. 
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Figure  8.25  shows  inverted  ratios  between  simulated  clear  atmosphere  and 
simulated  atmosphere  with  lOppmv  S02  divided  by  experimental  clear  atmosphere 
spectra.  There  is  a  significant  change  at  300nm.  Finally,  figure  8.26  illustrates  plots  of 
ratios  of  more  data  sets,  which  are  consistent  with  the  above  observations.  The  difference 
at  300nm  appearing  in  all  figures  is  caused  by  the  simulated  S02  concentration,  implying 
that  this  is  the  point  for  which  the  most  significant  S02  effect  in  measurements  should  be 
searched. 


RATIOS  OF  SELECTED  EXPERIMENTAL  TO  SIMULATED  fRES.OJnmt  DATA 
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Figure  8.24.  Ratios  of  data  taken  over  the  Sulfur  Works  vents  (SW_080  and  SW_206) 
and  clear  atmosphere  at  Lassen  Peak  (SW_503)  divided  by  simulated  clear  atmosphere. 
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RATIOS  OF  SELECTED  EXPERIMENTAL  TO  SIMULATED  (RES.  0.5nm)  DATA 
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Figure  8.25.  Inverted  ratios  of  simulated  clear  atmosphere  and  lOppmv  of  SO  2- 
containing  atmosphere  divided  by  clear  atmosphere  at  Lassen  Peak  (SW_503). 


RATIOS  OF  SELECTED  EXPERIMENTAL  TO  SIMULATED  (RES.0.5nmt  DATA 


Figure  8.26.  Ratios  of  selected  experimental  to  simulated  data. 
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Finally,  figure  8.27  compares  MODTRAN  with  S02  calibration  data.  The 
MODTRAN  curve  is  the  ratio  of  atmosphere  containing  lOppmv  SO  2,  divided  by  clear 
atmosphere.  The  S02  calibration  curve  is  the  ratio  of  the  data  acquired  with  lOOmmHg  of 
S02  in  the  photocell,  divided  by  the  data  with  the  photocell  evacuated.  There  is  a  very 
good  agreement  in  the  absorption  characteristics  of  the  gas  around  the  center  of  the 
maximum  spectral  response  of  the  filter.  Disagreement  away  from  300nm  is  due  to  the 
differences  between  real  and  simulated  filter  response.  There  is  also  a  shift  in  wavelength 
associated  with  the  focus  of  radiation  on  the  focal  plane  of  the  instrument  (Gray,  2002). 


characteristics. 


D.  SUMMARY 

Initially  MODTRAN  simulations  were  used  in  order  to  identify  the  absorption 
characteristics  of  S02  and  its  presence  signatures  at  various  concentrations  and  spectral 
resolutions.  Then,  experimental  data  processing  aimed  to  determine  if  those 
characteristics  and  signatures  were  present  in  the  images  taken  at  Sulfur  Works.  Finally,  a 
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combined  analysis  was  performed.  The  observations  deduced  from  the  analysis  are 
summarized  below: 

•  The  simulated  data  have  S02  absorption  characteristics  that  agree  very 
well  with  the  laboratory  absorption  cross  section  measurements  described 
in  chapter  III  (figures  3.2,  3.7  and  8.2  through  8.6). 

•  According  to  the  simulations,  S02  presence  in  the  atmosphere  causes 
differences  in  the  spectral  profiles,  which  are  very  small  in  small 
concentrations  but  increase  and  are  noticeable  in  larger  concentrations. 

•  Ratios  of  simulated  S02-containing  and  clear  atmosphere  reveal  a  clear 
and  characteristic  signature  even  for  low  concentrations  (figures  8,7  and 
8.8). 

•  The  spectral  profile  of  the  simulated  clear  atmosphere  in  the  region  of 
interest  is  mainly  due  to  solar  irradiance  (figure  8.12). 

•  According  to  the  S02  calibration,  LINUS  can  record  the  fine  spectral 
structure  caused  by  S02  absorption  bands  at  slit  widths  around  0.11mm 
(figure  8.14). 

•  Comparisons  of  the  experimental  data  taken  at  Sulfur  Works  (atmosphere 
possibly  containing  S02)  and  Lassen  Peak  (clear  atmosphere)  show  that 
the  most  noticeable  differences  occur  at  297.7  and  299.7nm  (figure  8.16). 

•  The  details  of  the  experimental  spectral  profiles  match  those  simulated  at 
0.5nm  resolution  much  better  that  those  of  the  finer  0.08nm  resolution. 
This  indicates  that  LINUS  with  a  0.11mm  slit  width  operates  at  about 
0.5nm  resolution  (figures  8.21  and  8.22). 

•  There  is  a  good  agreement  between  simulated  and  experimental  data  from 
about  290  to  306nm.  From  306nm  and  above,  the  experimental  profile 
becomes  wider  and  the  peaks  and  troughs  are  in  different  wavelength 
positions.  This  is  attributed  to  effects  of  the  filter  of  the  device  to  the 
recorded  radiation,  which  could  not  be  reproduced  by  the  filter  simulation. 

•  Combined  analysis  of  the  experimental  and  simulated  data  of  lOppmv  SO  2 
concentration  indicated  that  the  largest  noticeable  difference  is  at  299.9nm 
which  corresponds  to  the  maximum  of  the  absorption  cross  section  of  the 
S02  molecule  (figure  8.23). 

•  The  ensemble  analysis  approach  did  not  reveal  any  bold  S02  absoiption 
characteristics  in  the  experimental  data. 

The  final  conclusion  for  the  measurements  taken  at  Lassen  Park  is  that  there  are 
minor  indications  of  S02  presence  in  the  sulfuric  plums  emanating  at  Sulfur  Works.  The 
concentration  is  estimated  to  be  much  less  than  lOppmv.  This  result  is  in  contradiction 
with  what  was  expected  when  the  site  was  selected  for  the  LINUS  first  measurements. 
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Extensive  investigation  in  the  literature  in  order  to  find  out  any  previous  measurements  or 
gas  analysis  of  the  Sulfur  Works  emanations  did  not  reveal  anything.  The  U.S. 
Geological  Survey  was  conducted  and  communication  with  Mr.  M.  A.  Clynne  of  the 
Volcano  Hazards  Team  indicated  that  Sulfur  Works  plumes  should  not  have  any  or  have 
very  little  SCT.  This  is  the  case  for  any  gas  coming  through  an  aqueous  fluid,  because 
SO2  dissociates  into  HS  and  S,  and  in  effect,  is  scrubbed  from  the  plume.  This  is  in 
agreement  with  the  results  of  this  thesis. 
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IX.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  is  one  small  step  in  investigating  the  concept  of  passive  detection  of 
gases  in  an  extended  scene  in  the  atmosphere  by  using  the  principles  of  remote  s  ensing  in 
the  ultraviolet  region  of  the  electromagnetic  spectrum.  Specifically,  the  sensor  used  was 
LINUS  (Lineate  Imaging  Ultraviolet  Spectrometer)  and  the  gas  under  study  was  sulfur 
dioxide  (SCL).  Due  to  problems  with  the  scanning  software  at  the  time  this  thesis  was 
conducted,  LINUS  was  mainly  used  as  a  spectrometer  instead  of  an  imager. 

The  first  goal  of  this  project,  which  was  to  develop  the  operational  capability  of 
LINUS  and  deploy  it  for  the  first  time  in  the  field  in  order  to  acquire  data  of  natural 
plumes,  was  successfully  accomplished.  The  device  was  deployed  at  Lassen  Volcanic 
National  Park  in  August  and  September  of  2002,  and  a  number  of  measurements  of  the 
sulfuric  plumes  and  clear  atmosphere  in  the  area  were  acquired.  During  those 
deployments,  a  number  of  problems  were  encountered.  On  site  or  subsequent  solutions  to 
those  problems,  established  non-existent  knowledge  and  experience  in  using  the  device 
outside  the  laboratory.  The  response  of  the  device  was  also  investigated  by  changing 
various  settings,  such  as  the  slit  width,  acquisition  time  and  microchannel  plate  gain. 

Furthermore,  modeling  of  the  data  was  performed  using  MODTRAN4.  Modeling 
was  designed  to  simulate  the  signal  going  into  the  device  optics  section  after  it  had 
passed  through  the  UV  filter.  The  specific  conditions  under  which  several  measurement 
sets  were  taken,  were  simulated.  The  simulations  were  in  very  good  agreement  with 
laboratory  S02  absorption  cross  section  measurements  found  in  the  literature  and  gave  an 
indication  of  how  the  LINUS  data  should  look  like  in  the  wavelength-radiance  space. 
They  also  established  an  S02  absorption  signature  that  was  used  to  identify  S02 
absorption  effects  in  the  experimental  data. 

The  experimental  data  were  analyzed  using  an  individual  image  an  alysis  approach 
by  selecting  specific  images  and  investigating  their  spectral  structure,  and  an  ensemble 
approach,  where  the  images  were  processed  as  a  whole  in  order  to  identify  the  presence 
of  bold  SO2  characteristics.  The  first  approach  showed  minor  SCL  presence  indications, 
while  the  second  did  not. 
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Details  of  the  analysis,  findings  and  observations  are  presented  in  the  previous 
chapter.  The  final  conclusion  for  the  measurements  taken  at  Lassen  Park  is  that  there  are 
minor  indications  of  SO2  presence  in  the  sulfuric  plums  emanating  at  Sulfur  Works.  The 
concentration  is  estimated  to  be  much  less  than  lOppmv.  This  result  is  in  contradiction 
with  what  was  expected  when  the  site  was  selected  for  the  LINUS  first  measurements. 
Extensive  investigation  in  the  literature  did  not  reveal  any  previous  measurements  or  gas 
analysis  of  the  Sulfur  Works  emanations.  Communication  with  Mr.  M.  A.  Clynne  of  the 
U.S.  Geological  Survey  Volcano  Hazards  Team  indicated  that  there  should  be  little  or  no 
SO2  in  the  plumes  emanating  from  Sulfur  Works,  since  it  is  a  hot  springs  area  with  no  hot 
magma.  This  is  in  agreement  with  the  results  of  this  thesis. 

Follow-on  work  on  this  project  should  include  the  following: 

•  Replacement  of  the  existing  UV  filter,  with  a  new  filter  of  wider  spec  tral 
response  centered  at  300nm  and  higher  transmission  coefficient. 

•  Determination  of  LINUS  spectral  resolution  at  various  slit  widths  and 
measurement  of  the  SO2  detection  threshold  concentration. 

•  Selection  of  a  new  site  for  measurements  with  a  definite  and,  if  possible, 
known  or  estimated  S02  concentration  by  other  means,  so  that  direct 
comparisons  with  LINUS  results  can  be  done.  The  sites  could  be  either 
active  volcanoes  or  smoke  stacks  with  SO  2- rich  plumes. 
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APPENDIX  A.  MODTRAN  INPUT  CARDS  AND  VARIABLES 
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VARIABLE  VALUE  |  2.130  I  blank  |  70.000  I  0.000  I  0.000  I  0.000  I  0  I  0.000 
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APPENDIX  B.  EXAMPLES  OF  MODTRAN  INPUT  CARDS 


Two  examples  of  MODTRAN  input  cards  are  presented  in  this  appendix.  In  the 
first  three  pages,  the  input  data  are  presented  together  with  their  “coordinates”,  a  column 
number  from  01  to  80  in  the  first  and  last  two  rows  and  the  card  number  in  the  first  five 
columns.  These  “coordinates”  are  used  for  user/reader  convenience  and  do  not  exist  in  a 
normal  input  file.  In  the  last  two  pages  of  the  appendix,  a  normal  input  file  is  presented. 
The  files  simulate  measurements  at  Sulfur  Works  Location  2  with  a  lOppmv 
concentration  of  SO2  in  the  first  atmospheric  layer  with  thickness  90m. 


00000000011111111112222222222333333333344444444445555555555666666666677777777778 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

1  - M  732100000001  -1  0.000  0 

1A - TT  8F  0  365.00000  0  0  F  T  F  F  0.000 

1A2  —  DATA/BMP  9  9_01  .BIN 

2  -  100  00  0.00000  0.00000  0.00000  0.00000  2.13000 

2C -  36  1  0  PROTOTYPE  INPUT  CARD  -  10PPMV  OF  S02 

2C101  2.130  0.000E+0022.000E+00  0.000E+00  0.000E+00  0 . 000E+002B22222222A222  2 

2C201  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  1.000E+01  0.000E+00  0.000E+00 

2C201  0 . 000E+00 

2C102  2.220  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C202  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C202  0 . 000E+00 

2C103  2.976  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C203  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C203  0 . 000E+00 

2C104  3.732  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C204  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C204  0 . 000E+00 

2C105  4.488  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C205  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C205  0 . 000E+00 

2C106  5.244  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C206  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C206  0 . 000E+00 

2C107  6.000  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+00222222  22222222  2 

2C207  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C207  0 . 000E+00 

2C108  7.000  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C208  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000  E+00  0.000E+00  0.000E+00 

2C208  0 . 000E+00 

2C109  8.000  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C209  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C209  0 . 000E+00 

2C110  9.000  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 

2C210  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00 

2C210  0 . 000E+00 

2C111  10.000  0 . 000E+00  0.000E+00  0.000E+00  0.000E+00  0 . 000E+0022222222222222  2 
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2C211  0 . OOOE+OO  0 .  OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 
2C211  O.OOOE+OO 

2C112  11.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C212  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C212  0  .  OOOE+OO 

2C113  12.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C213  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C213  O.OOOE+OO 

2C114  13.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+ 0022222222222222  2 

2C214  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C214  O.OOOE+OO 

2C115  14.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C215  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C215  O.OOOE+OO 

2C116  15.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C216  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C216  O.OOOE+OO 

2C117  16.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C217  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C217  O.OOOE+OO 

2C118  17.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C218  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C218  O.OOOE+OO 

2C119  18.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C219  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C219  O.OOOE+OO 

2C120  19.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C220  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C220  O.OOOE+OO 

2C121  20.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C221  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C221  O.OOOE+OO 

2C122  21.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C222  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C222  O.OOOE+OO 

2C123  22.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C223  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C223  O.OOOE+OO 

2C124  23.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C224  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C224  O.OOOE+OO 

2C125  24.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C225  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C225  O.OOOE+OO 

2C126  25.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C226  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C226  O.OOOE+OO 

2C127  30.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C227  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C227  O.OOOE+OO 

2C128  35.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C228  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C228  O.OOOE+OO 

2C129  40.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C229  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C229  O.OOOE+OO 

2C130  45.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C230  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C230  O.OOOE+OO 

2C131  50.000  O.OOOE  +  OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C231  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

2C231  O.OOOE+OO 

2C132  55.000  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0 . 000E+0022222222222222  2 

2C232  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 
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APPENDIX  C.  SELECTED  SPECTRAL  PROFILES  FROM 
EXPERIMENTAL  DATA 


This  appendix  shows  plots  of  the  spectral  profiles  of  selected  images  taken  with 
LINUS  at  Sulfur  Works  Location  2  and  Lassen  Peak  (Location  3). 


SELECTED  DATA  FROM  SULFUR  WORKS  LOCATION  2  AND  LASSEN  PEAK 
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Figure  C.  1 .  Plot  of  the  spectral  profiles  of  all  selected  images. 
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SULFUR  WORKS  LOCATION  2  -  LOOKING  OVER  THE  VENTS 
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Figure  C.2.  Plot  of  spectral  profiles  of  one  group  of  selected  images  at  Sulfur  Works. 
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Figure  C.3.  Plot  of  spectral  profiles  of  another  group  of  selected  images  at  Sulfur 

Works. 
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Radiance  (Digital  Number)  Radiance  (Digital  Number] 


SELECTED  DATA  ACQUIRED  AT  LASSEN  PEAK  (LOCATION  3) 
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Figure  C.4.  Plot  of  spectral  profiles  of  selected  images  at  Lassen  Peak. 

DATA  ACQUIRED  WITH  COVERED  APERTURE  (NO  INPUT  SIGNAL1 


Figure  C.5.  Plot  of  spectral  profiles  of  measurements  taken  with  the  aperture  of 

LINUS  covered. 
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